Astronomy & Astrophysics manuscript no. 15338 


©ESO 2010 


September 9, 2010 





NGC6240: merger-induced star formation & gas dynamics* 

H. Engel 1 , R. I. Davies 1 , R. Genzel 1 , L. J. Tacconi 1 , E. K. S. Hicks 1 , E. Sturm 1 , T. Naab 2 , R H. Johansson 2 , S. J. Karl 2 , 

C. E. Max 3 , A. Medling 3 , and R R van der Werf 4 



1 Max Planck Institut fur extraterrestrische Physik, Postfach 1312, 85741 Garching, Germany 

2 Universitatssternwarte, Scheinerstrasse 1, 81679 Miinchen, Germany 

3 Center for Adaptive Optics, University of California, Santa Cruz, CA 95064, USA 

4 Leiden Observatory, Leiden University, P.O. Box 9513, NL-2300 RA Leiden, NL 



o 

<N 

Or 

00 



O 

U 

Or 

6 

•*-> 



> 
ON 

m 

i— H 

ON 

o 
o 



ABSTRACT 

We present spatially resolved integral field spectroscopic K-band data at a resolution of 0.1 3" (60 pc) and interferometric CO(2-l) 
line observations of the prototypical merging system NGC 6240. Despite the clear rotational signature, the stellar kinematics in the 
two nuclei are dominated by dispersion. We use Jeans modelling to derive the masses and the mass-to-light ratios of the nuclei. 
Combining the luminosities with the spatially resolved Bry equivalent width shows that only 1/3 of the K-band continuum from the 
nuclei is associated with the most recent star forming episode; and that less than 30% of the system's bolometric luminosity and only 
9% of its stellar mass is due to this starburst. The star formation properties, calculated from typical merger star formation histories, 
demonstrate the impact of different assumptions about the star formation history. The properties of the nuclei, and the existence of a 
prominent old stellar population, indicate that the nuclei are remnants of the progenitor galaxies' bulges. 

Key words, galaxies: active — galaxies: individual (NGC 6240) — galaxies: interactions — galaxies: kinematics and dynamics — 
galaxies: starburst — infrared: galaxies 



1. Introduction 

Major mergers are key drivers of galaxy evolution. Profoundly 
transformational events, they substantially affect virtually all 
properties of a galaxy. Mergers between disc galaxies are 
believe d to be responsibl e for t r iggering galaxy-wide star- 
bursts (iToomre & Toomrel Il972l iMihos & Hernquistl \\99t 



Genzel et all Il998a[lb 
etall 



2005; Hopkins 



, IVeilleux et all I20021 ISpringel elal 
2006I)), quasar activity (|S anders et al 1 
Sanders & Mirabell 19961: IVeilleux et all 120021: 
Springel et all 120051: iDasvra et all l2006bl 
and the creation of e llipti- 



1988a||b|; 
[ Jogeel 120041: 

Hop kins et all 120061" 120081). 
cal galaxies qToomre & Toomrel 119721: 
" Kormendv & Sandersl Il992l 



. Toomrel 
Mihos & Hernquistl 



1977 



1996; 



Genzel et all 120011 Das vra et alll2006allbl) . 

Star formation activity in merging gas-rich galaxies 
peaks between fi rst perip ass age and final co alescence 
(Sanders & Mirabel 1 19961 IMihos & Hernquistl 119961 
IVeilleux et all 120021 Springe l et all 120051) . Due to absorp- 
tion and re-emission by dust, almost all of this energy is 
emitted at infrared wavelengths, giving rise to extremely high 
infrared luminosities. Objects with Ltr = 10 11_1L9 L o are called 
Luminous Infrared Galaxies (LIRGs), those with Lir > 1O 12 L 
UltraLuminous InfraRed Galaxies (ULIRGs). The local LIRG 
population consists of bo th (typically early- stage) mergers and 
non-i nteracting galaxies (iSandersl Il992l lAlonso-Herrero et all 
120061) . whereas local ULIRGs are al most always i nteract- 
ing galaxies beyo nd first encounter (Sanders et all [l988a; 
San ders & Mirabell 119961 IVeilleux et all 120021: Ijogeel [2004). 
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* Based on observations at the European Southern Observatory VLT 
(079.B-0576). 



Local ULIRGs are predominantly powered by the sta rburst, 
except for the most luminous objects dGenzel et al.Lll998bh . 

Despite their recognised importance, we still lack a good un- 
derstanding of merger processes. This is due both to the com- 
plexity and the shortlived nature of these events. A rare oppor- 
tunity to study the transient phase between first encounter and 
final coalescence of two merging gas-rich spirals is afforded 
by NGC 6240, at a distan ce of 97Mpc. With Li R ~10 1L8 L o 
(ISanders & Mirabell 119961) . it falls just short of being formally 
classified as a ULIRG. HST images show a large-scale morphol- 
ogy dominat ed by tidal arms char acteristic for mergers past first 
peripassage dGerssen et all |2004). It is therefore more charac- 
teristic of the ULIRG class, to which it is commonly assigned; 
and our analysis shows its luminosity is likely to exceed the 
ULIRG threshold in the next 100-300 Myr. Two distinct nu- 
clei, with a projected separation of ~1.5" or 700 pc, are seen 
at infrared , optical, and r adio w avelengths (Tecza et al., 2000; 
Ma x et all 120051 120071 iGerssen et all l2004t iBeswick et all 
120011). Each nuc l eus i s host to an AGN, detected in hard X-rays 



2003 ) and at 5 GHz dGallimore & Beswickl 



Lutz et al. ( 2003J and|Xrmus et al. ( 2006) es 



( Kom ossa et all 
2004). However, 
timate that the AGN co ntribute less than half , and perhaps only 
25%, of the luminosity. iTacconi et al.l (Il999l) find the cold gas, 
as traced by CO(2-l) emission, to be concentrated between the 
two nuclei. The 1-0 S(l) H2 line emi ssion is one of th e most 
powerful found in any ga laxy to date ([J oseph et al., 1984 ), prob- 
ably excited in shocks ( van derWerfet all 119931: ISugai et all 



ll997HTecza etal.1.1 2000: 



Lutz et alll2003l) . The stellar kinemat- 



ics (iTecza et al. , 2000) display rotation centred on each of the 



two nuclei, and extraordinaril y large stella r velo city dispersions 
(-350 km s" 1 ) between them. ITecza et al] (2000) find evidence 
for a recent starburst i n the central kpc (the region encompass- 
ing the two nuclei), and Polla ck etall (12007) resolve a few dozen 
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clusters in the nuclear region that are also consistent with a re- 
cent starburst. 

In this paper we b uild on the qualitative picture put forward 
by iTecza et al.l (l200Qb that NGC 6240 is a merger already past 
its first close encounter, where the luminosity arises from a very 
recent burst of star formation, and therefore in which the mass 
must be due to an older stellar population. We use new adaptive 
optics integral field K-band spectra, interferometric mm CO(2-l) 
data, and merger star formation histories from numerical simula- 
tions, to perform a quantitative analysis of the mass, luminosity, 
age, and origin of the young and old stellar populations. 

We first introduce the observations and data reduction pro- 
cesses for the data sets used in our analyses (® . We briefly look 
at the molecular gas (® , and then focus on a detailed analysis 
of the stellar dynamics ^5] §0 before discussing the 
merger geometry and stage, and the curious CO(2-l) morphol- 
ogy (® . Bringing together these results, we then investigate the 
scale of the starburst (® and nature and origin of the nuclei 
( ClOl) . We summarise our results in CTT1 



2. Observations and Data Processing 




Fig. 1. K-band adaptive optics images of NGC 6240 in square- 
root colour scale. Left: K-band flux density from NIRC2 o n 
Keckll (pixel scale O.Or MPollack et al.ll2007l:lMax et al.ll2007l) . 
Right: 2.2 jum continuum from SINFONI on the VLT (pixel scale 
0.05"), covering approximately the same field as left. In both 
panels, circles indicate AGN positions (see north is up and 
east is left, and the red bar indicates 1 " (500 pc). 



2.1. SINFONI Data 



2.1 .1 . Observations and Reduction 

Observations of NGC 6240 were performed on the night 
of 20 Aug 2007 at the Very Large Telescope (VLT) with 
SIN FONI. SINFONI is a ne ar-infrared integral field spectrom- 
eter (lEisenhauer et all I2003I) which include s a curvature-based 
adaptive optics system (iBonnet et all I2004). For these observa- 
tions the wavef ront reference was provided by the Lase r Guid e 
Star Facility dBonacini Calia et all I20061 iRabien et all l2004l) . 
We used the R = 13.7 mag star 40" off-axis from NGC 6240 
to correct for tip-tilt motions. The instrument was rotated by 70° 
east of north in order to acquire this star. In order to cover both 
nuclei of NGC 6240 simultaneously, the 0.05"x0.10" pixel scale 
was chosen, giving a 3.2"x3.2" field of view. This pixel scale 
and the K-band grating provide a resolution of 6.6 x 10~ 4 //m, 
equivalent to R ~ 3300 or 90 km s" 1 FWHM at 2. 18 /mi. 

The data set comprised 12 object frames interleaved with 4 
sky frames. Each exposure was 300 seconds, yielding a total on- 
source integration time of 60 minutes. Data were reduced usin g 
the dedicated SPRED software package (lAbuter et all l2006h . 
which provides similar processing to that for long-slit data but 
with the added ability to reconstruct the data cube by stacking 
the individual slitlets. In order to optimise the sky subtraction, 
we made use of the algorithm described in IDaviesI (120071) . The 
image was rotated to restore the standard orientation, with north 
along the Y-axis, and east along the X-axis. The pixel scale of 
the final reduced cube was 0.05" x 0.05". 

Standard star frames were similarly reconstructed into cubes. 
Telluric correction and flux calibration were performed using 
HD 147550, a B9 V star with K=5.96mag. Velocity Standards 
were also observed with the same spatial and spectral settings, in 
order to enable recovery of the kinematics. The 3 stars observed, 
and added to our library of similar stars, were HD 176617 
(M3III), HD 185318 (K5III), and HD 168815 (K5II). These 
stars were chosen because of their deep CO absorption longward 
of 2.3 yum which provides a good match to that of NGC 6240. 



2.1.2. PSF Estimation 

We estimate the PSF by comparing the SI NFONI data to highe r 
resolution K-band Keck AO imaging data dPollack et all [2007). 
which are both shown in Fig.[T] Although the spatial resolution 
of our integral field observations is not as good as that achieved 
by these Keck imaging data, we note that the two bright spots 
in the northern nucleus are clearly separated in our data, indi- 
cating that our resolution is better than the 0.2" achieved by 
HST/NICMOS (Fig. 3c in lGerssen et aDl2004l) . 

As outlined by IDaviesI (I2008I) and employed by 
Mu eller Sanchez et~aU (I2006I) . since an observed image is the 
convolution of an intrinsic image with a PSF, I b s = I intr ®PSF, 
we can use a higher resolution image with a known PSF to 
obtain a good estimate of the PSF in a lower resolution image 
in two steps. Having resampled the data to the same pixel scale, 
we first find the broadening function B whose convolution 
with the high resolution image yields the best match to the 
low resolution image, Ii ow = high ® B. Since the underlying 
intrinsic image is the same for both l\ ow and high, we can then 
estimate the PSF of the low resolution image by convolving 
the broadening function B with the PSF of the high resolution 
image, PS F\ ow = PS Fhigh ® B. Since in this case B dominates 
the size and shape of our PSF, uncertainties in PSF high have 
little impact on PS F\ ow . For the analyses carried out here, this 
procedure gives a sufficiently accurate definition for the PSF. 
We find the PSF to be well matched by an asymmetric Gaussian 
with major axis oriented -20° east of north. The resolution 
is 0.097"x0.162" FWHM, corresponding to 50 x 80pc at 
the distance of NGC 6240. The asymmetry is to be expected 
since the resolution of the original data is pixel limited, being 
sampled at 0.05"x0.10". The reason we are able to do better 
than the formal Nyquist limit is because the SINFONI data are 
resampled to 0.05"x0.05"and many dithered exposures are 
combined using sub-pixel shifts to align them. 
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Fig. 2. The ratio between the measured dispersion <x ou t and the 
input dispersion cr in for a range of different stellar templates used 
to extract kinematics from a K 51b template. For each combina- 
tion, a variety of dispersions in the range 1 50-300 km s" 1 were 
used and noise was added to achieve a S/N-35. Variations in 
measured stellar kinematics using templates with a large range 
of Wco and stellar types are around a few percent, showing that 
with our method a mismatched stellar template has little effect. 



2.1.3. Spatial binning 

Because of the limited signal-to-noise in our data, we have 
binned them spatially using an optimal Voronoi tessellation 
( Cappellari et alll2003l) . This algorithm bins pixels together into 
groups by accreting new pixels to each group according to how 
close they are to the centroid of the current group. The result- 
ing groups then provide a set of positions (centroids) and mean 
fluxes which are used as the initial 'generators'. A further algo- 
rithm optimises the generator configuration based on a centroidal 
Voronoi tessellation. The final set of generators are the positions 
of the flux- weighted centroids of each binned group, and have 
the property that each pixel in the original image is assigned to 
the group corresponding to the nearest generator. This procedure 
only affects spaxels below a specified signal-to-noise threshold, 
and hence does not impact the spatial resolution of high signal- 
to-noise regions. The signal-to-noise cutoff was 20, chosen such 
that the regions with the highest signal-to-noise remained fully 
sampled, but lower signal-to-noise regions were binned. This 
avoided compromising the spatial resolution around the center 
of each nucleus, while enabling us to extend the region in which 
analyses can be performed. The effect of this can be seen in the 
resulting stellar kinematics maps, where the colours in the outer 
regions appear in blocks rather than individual pixels. The bin- 
ning scheme that this routine provided was applied to each spec- 
tral plane in the cube, and the kinematics were extracted from 
each bin as described below. 

2.1 .4. Extracting Stellar Kinematics 

The observed wavelength regime contains the stellar CO absorp- 
tion bandheads long ward of 2.29 yum, whose sharp blue edges are 
very sensitive to stellar motions. We utilise CO 2-0 and CO 3-1 to 
derive 2D maps of the stellar velocity and dispersion. The spec- 
tra were prepared by normalising them with respect to a linear fit 
to the line-free continuum. The continuum level was then set to 
zero by subtracting unity. A template spectrum, that of the M3 III 
star HP 176617, w as p repared in the sam e way. Although both 
ISugai et al.1 (1 1997b and lTecza etafl (l2000h used KIb stars, they 
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Fig. 3. Spectra (black lines) and template fits (green lines) to the 
CO 2-0 and CO 3-1 stellar absorption features. These are for in- 
dividual bins in the southern nucleus (top), the northern nucleus 
(middle) and in the region of high dispersion between them (bot- 
tom). The blue line indicates the continuum (i.e. zero) level; red 
asterisks denote pixels rejected during the fit; the pairs of ma- 
genta lines enclose regions around the bandheads for which the 
weighting is enhanced during the fit. These plots demonstrate 
the quality of the fit using our single M3 III template. 



also both noted that Mill stars provide d almost equally good 
fits to the K-band features (see Fig. 5 in lTecza et alj|2000b . We 
therefore do not expect template mismatch to be a problem but, 
as explained below, our method furthermor e minimises the im- 
pact of any discrepancy. ISilge & Gebhardtl (120031) showed that 
the dispersion measured depends on the equivalent width of the 
template; this can be understood since the depth of the extinction 
feature relative to the continuum level (which may contain non- 
stellar emission) places a strong constraint on the fitting param- 
eters. But it is the spectral width, not the equivalent width (i.e. 
depth), of the absorption feature that is important for deriving the 
dispersion. In our fitting method we resolve the problem of the 
feature's depth by setting the continuum level to zero. Therefore 
matching the relative depth of the absorption features, and hence 
the choice of template, is no longer critical. 

In order to demonstrate this q uantitatively, we ha ve made 
use of the stellar template library of lWinge et al] (l2009h . We cre- 
ated a large series of mock observations by broadening a K5 lb 
template by various amounts in the range 150-300 km s" 1 , and 
adding noise. We then extracted the dispersion using various 
mis-matched stellar templates spanning a wide range of Wco 
(equivalent width of the CO absorption feature) and spectral 
type. As can be seen in Fig. [2J our method recovers the input 
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Fig. 4. Stellar velocity and corresponding noise map (top and bottom left, both in kms with respect to systemic), dispersion and 
corresponding noise map (middle, both in kms -1 ), and Wco 2-0 and its noise map (right, both in A). Contours trace the continuum 
as in Fig.[T] Black/magenta circles mark the AGN positions (see §0). The bar indicates 1" (images are 2.45 // x3.95 // ), the PSF size 
is shown in the lower left corner of panel 1 . North is up and east is to the left. 



dispersion with a typical error of only a few percent, regardless 
of the template's spectral type. This shows that our method of 
extracting the stellar kinematics works reliably even if the tem- 
plat e used is not a p erfe ct fit. We emphasi ze that, in contrast 
to lSugaiet all (1 19971) and lTecza etaD d2QQQh . we do not use the 
quality of the template fit to draw conclusions about the preva- 
lent stellar population. This requires caution, and we discuss the 
pitfalls associated with doing so in §0 As such, because we use 
this star only for extracting kinematics, our choice of template 
can be fully justified by the good fits to the bandheads in the 
galaxy spectra in Fig. [3] 

For the spectrum at each spatial pixel ('spaxel' hereafter) in 
the data cube, we convolve the template with a Gaussian, adjust- 
ing the properties of the Gaussian to minimise the difference be- 
tween the convolution product and the galaxy spectrum. Deviant 
pixels are rejected from the fit. The spectral regions covering 
the steep edges of the bandhead at all expected stellar velocities 
are given five-fold weight (Figure [3]), in order to focus the fit 
on the kinematics rather than details in the spectrum. To ensure 
consistency in our analyse s, the spectral ranges app lied to com- 
pute Wco were taken from lForster Schreibeil (l2000l) . since these 
ranges are also used by the stellar synthesis co de STARS which 
we u se later. A high signal-to-noise (> 50, Cappe llari et al.l 
2009) is required in order to measure deviations in the line of 
sight velocity distribution from a simple Gaussian. Without bin- 
ning beyond a useful spatial resolution, the signal-to-noise in 
most of our spectra is not sufficient to include the Gauss-Hermite 
terms I13 and in the fit. We therefore limited the kinematic 
extraction to V and cr. The noise in the kinematics was deter- 



mined using Monte Carlo techniques: 100 realisations of the 
data were generated, a number that ensures a fractional uncer- 
tainty in the noise of less than 10%, by perturbing the flux at each 
pixel according to its RMS. These were fit using the same pro- 
cedure as above. Although the RMS at each pixel can be derived 
when combining the data cubes, calculating the noise in bins is 
extremely difficult due to the correlatio ns between neighbour- 
ing pix els. This is discussed in detail by iForster Schreiber et al.l 
( 2009). Instead, we have derived the noise as the RMS of the dif- 
ference between a spectrum and the convolved template (exclud- 
ing rejected pixels). This may slightly over-estimate the noise, 
but will then yield conservative error estimates on the kinematic 
parameters. The resulting final maps of stellar velocity, disper- 
sion, and Wco 2-0 are displayed in Fig. HI 

2.2. Plateau de Bure Interferometer Data 

We have mapped the 12 CO(J = 2-l) line in NGC 6240 with the 
IRAM millimetre interferometer, which is l ocated at an alti- 
tude of 2550 m on the Plateau de Bure, France ( Guilloteau et al., 
Il992l) . The data were obtained in May 2007. The array consisted 
of six 15 m antennae positioned in two configurations providing 
a large number of baselines ranging from 32 to 760 m. We ob- 
served NGC 6240 for 5 hrs in A configuration, which provides 
a nominal spatial resolution of 0.35" at 230 GHz. However, the 
low (+2deg) declination of NGC 6240 limited the resolution in 
an approximately north- south direction. A spectral resolution 
of 2.5 MHz, corresponding to 3.3 kms -1 for the CO(2-l) line, 
was provided by 8 correlator spectrometers covering the total re- 
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Fig. 5. Map of CO(2-l) flux with spectral line profiles shown for selected spaxels. Circles denote AGN positions, stars denote 
positions where spectra were extracted. The beam size is displayed in lower left. The bar denotes 1", and north is up and east is to 
the left. 




Fig. 6. Map of H2 1-0 S(l) emission, with line shapes at selected spaxels. Circles mark AGN positions, stars denote positions where 
spectra were extracted, white contours represent continuum emission. The bar denotes 1", and north is up and east is to the left. 



ceiver bandwidth of 1000 MHz (1300km s" 1 ). All the data were 
first calibrated using the IRAM CLIC software. We then made 
uniformly weighted channel maps for the CO(2-l) data. We 
CLEANed all the maps using software available as part of the 
GILDAS package. To increase the sensitivity, maps were made 
with a velocity resolution of 26.6 km s" 1 . The CLEANed maps 
were reconvolved with a 0.69 // x0.26 // FWHM Gaussian beam. 
The rms noise after CLEANing is 1.4 mJy beam -1 . Fig.[5]shows 



the velocity-integrated CO(2-l) emission (black circles denot- 
ing approximate AGN locations), velocity moment map, and the 
dispersion (obtained from fitting Gaussians to the line shape). 
The AGN positions were es timated using the positions given by 
Gallimore & Beswick (2004); due to the astrometric uncertain- 
ties, these are only accurate to within 1-2 pixels (0.7-0. 14"). 
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n 



Fig. 7. Velocity moment (left) and dispersion (right) maps of H2. 
Magenta circles mark AGN positions, white contours represent 
H2 1-0 S(l) line emission. 



3. Molecular Gas Emission 



3.1. H 2 1-0S(1) 



NGC 6240 has 
any galaxy to 



the s trongest 
date 



Ho 



l ine e mission found in 
dJoseph et all Il984). Several authors 
(Ivan der Werf et allll993l:ISugai et allll997tlTecza et all I2Q0C 



)0) 
al. 



identify shocks as the excitation mechanism, and lOhyama et 
( 2003) use the relative intensity of H2 and CO(2-l) emission to 
conjecture that the shocks occur due to a superwind outflowing 
from the southern nucleus colliding with the molecular gas con- 
centration. Our new SINFONI data provide the most detailed 
spatial and spectral view currently available. As can be seen in 
Fig-El the line profiles are highly complex, with clear evidence 
of multiple components. This, and the dispersion map (Fig. [7]), 
support the picture of highly disturbed, turbulent gas. The emis- 
sion morphology is uncorrelated with the stellar emission dis- 
tribution, but roughly consistent wi th th e cold gas emission (as 
traced by the CO(2-l) emission, S3. 2b . The velocity moment 
map (Fig. [7]) shows a global velocity gradient, again broadly con- 
sistent with what is seen in CO(2-l). 

3.2. CO(2-1) 

Fig.[5]shows the velocity-integrated CO(2-l) emission, with line 
shapes at selected positions. With a factor of two improvement 
in reso lution along the E-W direction, our data re-affirm the find- 
ings of iTacconi et al.l (Il999h . However, our resolution in the cru- 
cial N-S direction is limited by the declination of the source. As 
a result we are not able to draw more detailed conclusions about 
the CO(2-l). 

The CO emission is concentrated in between the two nuclei, 
and the CO(2-l) velocity map, red and blue wings, and P-V- 
diagram all display sig ns of a velocity grad ient. This is consis- 
tent with the findings of lTacconi et al.l (Il999l) . who conclude that 
the molecular gas likely is concentrated in a self-gravitating ro- 
tationally supported, but highly turbulent disc in the internuclear 
region. However, such a central gas concentration is not expected 
from merger simulations, which generally predict the gas to re- 
main largely bound to the progenitors until the nuclei coalesce. 
We discuss this issue further in §[8] 

The total flux i s about 40% less than that measured by 
Tacconietal. (1999), indicating that our smaller beam has re- 
solved out some emission on intermediate to large scales. We 
have not made a correction for this, since we are interested pri- 



wavelength [um] 

Fig. 8. Comparison of different extinction models: We measured 
HST photometric data points (red squares) in 0.2" radius aper- 
tures; positions for the three examples shown here are indicated 
in the left panel. To these we fitted a synthetic STARS spec- 
trum which was reddened using different reddening prescrip- 
tions: screen extin ction (black), mixed extinction (grey), and the 
ICalzetti et al.l (l200Qh reddening law (magenta). We also overplot 
our SINFONI spectrum (blue), scaled to the HST 2.22mum data 
point (absolute SINFONI flu x at 2.2mum show n as blue dia- 
mond). As can be seen, the lCalzetti et al.l (l2000h reddening law 
provides the best fit to the data. 



marily in the small scale emission of the nuclei themselves, 
where the effect is likely to be negligible. To estimate the gas 
mass associated with the nuclei, we have measured velocity- 
integrated line fluxes within the same radii as those used in our 
dynamical mass modelling (§|7]), which were chosen to cover the 
extent of observed stellar rotation in the nuclei. In the northern 
nucleus we find S coA V = 37 Jy km s" 1 out to a radius of 250 pc; 
in the southern nucleus we find Sco^V = 178 Jy km s" 1 out to a 
radius of 320 pc. 

We estimate the gas mass by converting to line lu- 
minosity L'co = 3.25 x 10 7 S C oAVv-^d£(1 + z )-3 ( with 
£> L =97Mpc, z = 0.0243, v ob s,co 2 , = 225. 1 GHz), and using 



the relations L(CO 2 -i)/L(COi_ ) = 0.8 dCasoli et all fl9 92) 
and Mfa/Lco^ ~ 1 M /(Kkm s" 1 pc 2 ), as iDownes & Solomon 
(119981) find for ULIRGs. This yields ~ 0.2 x 10 9 M o and ~ 
1.1 x 10 9 M o for northern and southern nuclei, respectively. 
Similarly, for an aperture of 1" diameter centred on the CO emis- 
sion peak, we derive a gas mass of ~3.1xl0 9 M o . We caution 
that the nuclear masses are very uncertain and should be treated 
only as order-of-magnitude estimates. The reason is that a sig- 
nificant fraction of the gas within the apertures may not be phys- 
ically co-located with the nuclei, which would imply the masses 
are upper limits. On the other hand, the CO abundance around 
the nuclei is very likely to have been significantly reduced by X- 
ray irradiation from the AGN (see ® . In this case, the masses 
would be underestimates. 



4. Extinction and Luminosity of Nuclei 

Like most ULIRGs, N GC 6240 contains significant amounts of 
dust (iTecza et alll2000h . and hence correcting any flux measure- 
ment for extinction is paramount to ensure accuracy in the anal- 
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yses described below. It is valuable to consider the extinction de- 
rived from mid-infrared data to assess whether the near infrared 
might be affected by saturation. However, this appears to be un- 
certa in - primarily be cause no H n lines were detecte d either by 
ISO dLutz et all [20031) or Spi tzer dArmus et al.ll2006l) . Based on 
silicate absorption at 9.7 /mi. lArmus et al.l (120061) find an extinc- 
tion of Ay ~ 95 mag to the c oronal line region . The more mod- 
erate estimates presented by Lutz etaD (2003) are for a global 
extinction, and correspond to a dust screen Ay -15-20 mag. 
This estimate is comparable to der ivations at radio (Nh ~(1.5- 
2)xl0 22 cnr 2 ,lBeswick et~ai1 l200ll) and X-ray (N H ~ 10 22 cm" 2 , 
iKomossa et al. 120031) wavelengths. This is important, because 
it implies a modest typical K-band screen extinction of only 
A K ~ 1-2 mag. Thus our data should be sensitive to the majority 
of the emission. However, we note that the X-ray spectrum at en- 
ergies above lOkeV indicates that at least one of t he AGN them- 
selve s may be obscured by a column > 1 24 cm " 2 dVignati et all 
119991: llkebe et all l2000l: iNetzer et all 12005 ). Differences be- 
tween these various column density measurements are to be ex- 
pected since they sample different spatial scales and sight lines; 
for example, the X-ray data specifically measure gas column 
density along the line of sight to the AGN rather than the ex- 
tended cold dusty medium. 

When correcting for the obscuring effect of dust, one needs 
to make assumptions about its location and distribution. The 
two most commonly used models assume either a uniform dust 
screen between the observer and the stars, leading to a reduc- 



tion in observed flux according to F b s /F ei 



= e 



Or the 



stars and dust are assumed to be spatially coincident and uni- 
formly mixed ('mixed model'), in this case the observed flux 
decreases as Fob s /F e m = (1 - e~ T )/r. Another alternative is the 
ICalzetti et al.l (l2000l) reddening law, derived empirically from 
observations of starburst galaxies. Its wavelength dependence 
reflects both the dust grain properties and the distribution of 
the dust with respect to the stars in these galaxies. The effect 
of this reddening law is remarkably similar to a combination of 
mixed and screen extinction, with - for the degree of extinc- 
tion seen in NGC 6240 - the mixed component dominating in 
the infrared, and the screen component increasingly important 
at optical wavelengths. In order to investigate which extinction 
model best captures the characteristics of the dust distribution 
of NGC 6240, we obtained a number of archival HST imaging 
data, spanning wavelengths from 0.45/mi to 2.22/mi, and ex- 
tracted photometric data points from a number of 0.2" radius 
apertures across our field of view. And we calculated a syn- 
thetic stellar spectrum using the stellar synthesis code STARS 
( Sternberg, 1998; Sternberg et al., 2003; Forst er Schreiber et all 
120031: iDavies et all I2003L I2005L I2006L 12007b . assuming a star 
formation rate typical for a merger (Antennae' -simulation, ^8]). 
We then fitted this spectrum to the HST data points and our 
SINFONI K-band s pectra, using t he scr een and mixed extinc- 
tion models, and the Calze ttiet all (12 000) reddening law. As can 
be seen in Fig. [U the ICalzetti et all (l200Qh reddening law best 
reproduces the observations; the mixed model saturates at opti- 
cal wavelengths and cannot redden the spectra sufficiently, and 
the screen model is reddening the spectra too much at shorter 
wavelengths. We also test whether our choice of star formation 
history has an influence on this result, by conducting the same 
test with synthesised spectra for a 20 Myr old instantaneous star- 
burst and 1 Gyr o f continuous star formation. Both also indicate 
the ICalzetti et al.l (|2000) reddening law to be appropriate. 

We can then find F b s /F em for each spatial pixel by adjust- 
ing the redde ning for a set of st ellar template spectra according 
to Eqn. 2 in ICalzetti et al.l (l2000l) . until the best-fit of a linear 



combination of templates to the measured line-free continuum 
is achieved. Fig. [9] shows the resulting map of A K . The choice of 
stellar templates does not affect the result, since the K-band sam- 
ples the Rayleigh- Jeans tail of the blackbody curve and hence to 
a good approximation all stars, and late-type stars in particular, 
have the same spectral slope. 

One might ask whether a component of non-stellar emission 
from hot dust might be present, in our field of view in general 
and close to the AGN in particular. This is an important ques- 
tion, since it would have an influence on the stellar masses and 
mass-to-light ratios we calculate later on. That the K-band con- 
tin uum may even be d ominated by hot dust has been proposed 
by lArmus et all (120061) . However, this conclusion was based on 
a fit to the near-IR spectral energy distribution that was con- 
strained primarily as the residual in the blue side of a much more 
dominant cooler component, and is therefore rather uncertain. 
Instead, the good match of reddened synthetic stellar spectra 
to the photometry (Fig. [8]) indicates that non-stellar emission is 
unlikely to contribute significantly. More importantly, emission 
from hot dust would dilute the stellar CO absorption features. 
Since we measure Wco 2-0 ~ 12- 13 A (Fig.©, and the theoreti- 
cally possible (achievable only through a 10 Myr old instanta- 
neous starburst) maximum is ~ 1 8 A, this puts a firm theoretical 
upper limit of <30% on any non-stellar contribution, and realis- 
tically makes anything larger than a few percent unlikely. And 
since we do not see a localised dip in Wco 2-0 around the AGN 
positions, we conclude that any hot dust emission associated 
with the AGN is completely obscured at near-IR wavelengths. 

Table Q] lists the observed and dereddened luminosities mea- 
sured within apertures of diamet er \" cen t ered o n the north- 
ern and southern nuclei, from iTecza et all (|2000), this work, 
and from archival NICMOS data; and also those obtained by 
integrating the luminosity profiles out to 250 pc. We included 
the NICMOS data to obtain a third independent measurement, 
si nce our measurem ents of the observed luminosities and those 
of ITecza et all (l200Qh differ by a factor of two. We find the 
NICMOS data agree with our mea surements to within 10%. As 
can be seen, for lTecza et all {2000), the dereddening only alters 
the measured luminosities by -10%, whereas for us, the differ- 
ence is circa a factor of two. This difference is most likely due to 
the fact that we applied a spatially dependent, rather than single- 
valued, correction. 



5. Kinematic Centres and Black Hole Locations 

In this section, we attempt to confirm the hypothesis that the 
black hole locations do identify the centres of the progenitors, by 
independently determining the locations of the kinematic centres 
from the stellar velocity field. As we outline below, for a num- 
ber of reasons locating the kinematic centres of the observed 
stellar rotation reliably and accurately is extremely difficult to 
do. Neverth eless, within the uncertainties, we find that the BH 
positions of iMax et aH (120071) are consistent with the kinematic 

ce ntres. 

IMax et al.l (Hp07) determined the positions of the two AGN 
in NGC 6240 by combining images taken in the near-infrared, 
radio, and X-ray regimes. Postulating that the southern sub-peak 
of the n orthern nucleus as se en in the K-band ('NT in the no- 
tation of iGerssen et al.ll2004l) is coincident with the position of 
the northern black hole, they superpose radio data to infer that 
the southern AGN is located to the north-west of the 2.2 /mi 
peak of the southern nucleus. The offset is explained in terms 
of dust obscuration, which is supported by their 3.6//m images 
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Table 1. Integrated Luminosities of the Nuclei 
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North 


South 
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dereddened 
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[10 8 L o ] 


[10 8 L o ] 
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Apertures, Tecza et al. (2000) 


2.2 


2.3 


5.8 


6.8 


Apertures, this work 


5.6 


9.8 


15 


40 


Integrated profiles, this work 




11 




59 


Apertures, NICMOS 


5.3 




17 





References. Apertures are all diameter 1", the profiles are integrated out to a radius of 0.5". 




Fig. 9. Effective K-band extinction A#. As for Fig. \T\ contours 
trace the continuum and black circles denote the AGN positions. 
Each 0.05" pixel length corresponds to 25 pc at the distance of 
NGC 6240. 



showing the southern continuum peak to be coincident with the 
posited black hole location. Our data also support the existence 
of higher obscuration at this location: Fig. [9] shows that the re- 
gion suffering the greatest extinction overlaps with, and extends 
to the north west of, the south ern nucleus. We identify the AGN 
positions of M ax et all (|2007) to an accuracy of better than ± 1 
pixel (0.05") on our data using the 2.12/mi continuum features 
as reference points (as shown in Fig.Q]), and the angular separa- 
tion of the VLBA radio sources ( 1.5 11 +0.003". iMax et alJ l2007l 
and references ther ein). The uncertai nties on the BH positions 
as identified on the IMax et all (120071) data are smaller than the 
uncertainty introduced by this translation ont o our data. 

We use kinemetry ( Kra jnovic et all 12006) to parametrize the 
stellar velocity field which, as Fig. [4] shows, clearly exhibits or- 
dered rotation around each of the two nuclei separately. In many 
cases, the centre of a galaxy and its axis ratio (or inclination) 
and position angle (PA) can be extracted straightforwardly from 
the isophotes. This is not possible for NGC 6240 because the 
K-band isophotes are strongly asymme t ric w ith respect to the 
black hole locations found by lMax et al.l (120071) . Hence we must 
use our analysis of the velocity field to also fix these parameters 
simultaneously. 

Kinemetry decomposes a velocity map K(a, if/) into a series 
of elliptical rings that can be expressed as the sum of a finite 
number of terms: 

N 

K(a, i/O = Ao(a) + ^ A n (a) sm(ni//) + B n (a) cos(nif/). 

n=\ 



By minimizing particular terms, one can in principle find the 
centre, position angle, and axis ratio of the velocity field. Under 
the assumption that the velocity field is due to an axisymmetric 
thin disc, these correspond directly to the equivalent parameters 
for the galaxy (although one should bear in mind that, for exam- 
ple, with non-circular orbits the kinematic and isophotal major 
axes may not coincide; and for a geometrically thick system the 
axis ratio may not represent the inclination). 

The correct choice of kinematic centre minimises A2 and B2 
(which are also weakly dependent o n ellipticity and PA) and to 
a lesser degree Ai, A3, and #3 (see iKrajnovic et all 120061 for a 
more intuitive description of these parameters). The correct PA 
yields minimal values of A 1, A3, and By, correct ellipticity min- 
imises B3. For a range of different ellipticities and PA, we place 
the kinematic centre at each pixel within a grid of 11x11 pixel 
centered on the AGN location and calculate the corresponding 
sum of A2 and B^. This results in a set of 2D maps of A2 + #2 for 
each combination of PA and ellipticity. The position of the min- 
imum of each map then yields the best-fitting kinematic centre 
for that specific combination of PA and ellipticity. If the best-fit 
to the kinematic centre is the same for different values of ellip- 
ticity and PA (i.e. a unique best estimate for the kinematic centre 
exists), we then proceed to determine A 1 and A3 at this kinematic 
centre for a range of different position angles, taking as the best 
estimate of the PA that value which results in A\ + A3 being min- 
imal. Finally, the ellipticity is found by analogously minimising 
the B3 coefficient. 

With this method we are able to find a unique solution for 
all three parameters for the southern nucleus. But it fails for the 
northern nucleus where the best estimates for kinematic centre 
and position angle are interdependent: the minima for A2 + #2 
(determining the kinematic centre) lie on an arc with the exact 
position of the kinematic centre dependent on the choice of po- 
sition angle. The only way out of this impasse is to make ad- 
ditional assumptions, and so we arrive at a 'best-fit' by assum- 
ing the velocity field to have the same PA as the major axis of 
the continuum emission. We note that the PA of the continuum 
emission is independent of the extinction correction, because the 
extinction and continuum maps have the same PA. Figures \W\ 
and [TT] show the observed velocity field, the model, and the 
residuals for 2 cases in each nucleus: when the best- fitting centre 
derived from the kinemetry is used, and when the centre is fixed 
at the location of the black hole. 

In the northern nucleus, the location of the derived kinematic 
centre is reasonably consistent with the black hole location, dif- 
fering by only 0.12". Since the average residuals per pixel are 
33.5 km s" 1 and 35.3 km s" 1 for these two cases respectively, 
they can be considered statistically indistinguishable. We con- 
clude that in the northern nucleus, we can confirm that the black 
hole is located at the kinematic centre of the stellar rotation. 

For the southern nucleus, the difference is more significant. 
The separation between the best-fitting centre and the AGN lo- 
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Fig. 10. Northern Nucleus: Comparison of kinemetry models for 
2 different PAs and kinematic centres. It should be noted that the 
rotation curves do not differ significantly for the two different po- 
sitions of the kinematic centre. Black circles indicate black hole 
positions; black stars represent the kinematic centre; the dashed 
white line traces the major axis of rotation. Left to right: veloc- 
ity map, model, and residuals. First row: best-fitting kinematic 
centre from kinemetry analysis. The average residuals per pixel 
are 33.5 km s" 1 . Second row: kinematic centre fixed at the posi- 
tion of the AGN. The average residuals per pixel are 35.3 km s" 1 . 
Image size is 820x720 pc. 



cation is 0.22". The mean residuals are 28.2 km s" 1 (best- fit) and 
45.9 km s" 1 (AGN position). This can be understood with refer- 
ence to the dispersion map in Fig. [4] and the extinction map in 
Fig. [9j The region to the northwest of the AGN in the southern 
nucleus exhibits both anomalously high dispersion and extinc- 
tion - features that are to be expected in a merger system. It is 
therefore not clear whether the velocity field in this region is 
really tracing the rotation of the progenitor, or something more 
complex. Inspection of the residuals in Fig. Qj] suggests the lat- 
ter, an issue that we discuss in more detail in §0 Our conclusion 
for the southern nucleus is that the centre we derive from the 
velocity field is biased by perturbed kinematics on the northern 
side of the nucleus. 

In both nuclei, for further analyses we therefore adopt the 
black hole locations of iMax et~aD (120071) as identifying the kine- 
matic centres of the nuclei, and make use of the rotation curves, 
dispersion, and luminosity profiles centered at these positions. 
We derive axis ratios of 0.70 ± 0.02 and 0.65 ± 0.02 and PAs of 
229°±2° and329°±2° (measured east of north, PA pointing from 
receding to approaching velocities) for the northern and south- 
ern nuclei, respectively. We interpret the axis ratios in terms of 
an inclination for a flat system, noting that this may be an over- 
estimate since the nuclei are likely to be thick. The impact is that 
we may also overestimate the intrinsic velocity. 



6. Stellar Kinematics between the Nuclei 

NGC 6240 is known to exhibit an exce ptionally large stellar ve- 
locity dispersion between the nuclei (iLester & Gaffneyl Il994t 
iDovon et all 1 19941: iTecza et all l2000h . In Fig. g] we present a 
detailed 2D dispersion map which shows that the region with 
the highest (>300kms _1 ) dispersion is fairly localised, and lies 
across the northern side of the southern nucleus. The localised 
nature of this reg ion, together with p oor signal-to-noise, most 
likely explain whv lGenzel et al.l (l200ll) and lDasvra et al.l (l2006bl) 



1 1 1 B 1 1 1 1 1 1 
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Fig. 11. Southern Nucleus: Comparison of kinemetry models for 
different PAs and kinematic centres. It should be noted that the 
rotation curves do not differ significantly for the two different po- 
sitions of the kinematic centre. Black circles indicate black hole 
positions; black stars represent the kinematic centre; the dashed 
white line traces the major axis of rotation. Left to right: veloc- 
ity map, model, and residuals. First row: best-fitting kinematic 
centre from kinemetry analysis. Average residuals per pixel are 
28.2 km s" 1 . Second row: kinematic centre fixed at the position 
of the AGN; average residuals per pixel are 45.9 km s" 1 . Third 
row: kinematic centre fixed at position of the AGN and per- 
forming fit only to southern half of nucleus. Average residu- 
als per pixel over southern half are 21.1 km s" 1 . Image size is 
600x840 pc. 



cite a significantly lower maximal dispersion, as it can easily be 
missed by slit measurements. 

This high dispersion cannot be associated with the estab- 
lished stellar population in the southern nucleus because the high 
velocity (-500 km s" 1 at 500 pc) implies a dynamical timescale 
of tdyn ~5 Myr. Any asymmetry in the dispersion of a stellar pop- 
ulation orbiting the southern nucleus would be dispersed within 
this timescale. An alternative explanation is that this represents 
a region where emission from stars that have been formed re- 
cently as a result of the interaction is superimposed on the light 
from the nucleus. As the stellar kinematics are derived from ab- 
sorption lines, a superposition of two populations at different 
line-of- sight velocities would lead to an overestimate of the dis- 
persion: the line resulting from such a superposition would ap- 
pear broader than the intrinsic line widths of each of its con- 
stituent lines. Fitting the resulting line profile with a single kine- 
matic component would lead to a significant over-estimate of the 
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dispe rion. Support for this explanation is lent by lQhyama et al.l 
(2003). Based on a study of the morphology and kinematics of 
the warm and cold gas, they argued that clouds along the line 
of sight to the northern half of the southern nucleus were being 
crushed in the interaction. In this scenario, if the cloud crushing 
leads to star formation, there would indeed be two superimposed 
populations at this location. This is supported by Spitzer obser- 
vations of the Antennae galaxies, which show that the largest en- 
ergy output, and thus a high rate of star formation, occurs in the 
overlap region between the two galaxies' nuclei (Br andl etall 
[2QQ9h . 

Although rather speculative, in the following discussion we 
consider whether this scenario can in principle account for both 
the increased dispersion and irregularities in the velocity in the 
northern half of the southern nucleus. Our aim is to use a very 
simple toy model to test the basic validity of the hypothesis by 
reproducing the characteristic features of the data; we do not 
attempt to match them exactly, nor to constrain the numerous 
parameters that would be required to do so. 

We construct a two-population model as follows: the disc 
population is represented by a disc with the kinematic centre at 
the position of the AGN. PA, inclination, and velocity curve are 
found via a minimisation such as best to match the observed ve- 
locity field in the southern half. The dispersion is adopted from 
the radial dispersion curve derived from the southern half. We 
then subtract the disc velocity field from the observed veloc- 
ity field in the northern half, and take this to be the velocity 
field of the second stellar population (thus ensuring that the ob- 
served velocity field is reproduced as closely as possible). In or- 
der to derive the resulting dispersion, we combine the two popu- 
lations at each pixel by adding two Gaussians, with their centres 
corresponding to velocity and FWHM corresponding to disper- 
sion (resembling the superposition of two absorption or emis- 
sion lines). The resulting velocity field of the combination of 
both populations is then extracted by fitting a Gaussian to the 
result (resembling the method through which stellar kinematics 
are derived). As can be seen in Fig. [T21 the model, whilst not 
matching the observations exactly, does reproduce the charac- 
teristic features of the data, displaying locally significantly in- 
creased velocity dispersions. The fact that the dispersion map 
is not reproduced exactly is due to the fact that we intention- 
ally kept the model simple, assuming spatially constant, equal 
weightings and line FWHMs. We emphasise that the purpose of 
this exercise is simply to show that it is in principle possible to 
account for the locally significantly increased inferred dispersion 
through the effect of fitting a single velocity and dispersion to an 
absorption feature to which two stellar populations at different 
relative velocities have contributed. Although our assumptions 
that the two populations have the same luminosity and disper- 
sion are simplistic, this underlines the point that a significantly 
better fit could be achieved if these parameters were allowed 
to vary. However, we feel that pursuing this is unjustified given 
constraints available from the data. 

The total extinction-corrected luminosity emitted by the in- 
ternuclear area with cr > 330km s" 1 is 8.5 x 1O 8 L . If this can 
be attributed in roughly equal shares to the old stellar popula- 
tion and the newly formed stars, it would imply a stellar mass of 
~ 1O 8 M contained in the starburst population in this area. We 
speculate that this population may have originated from star for- 
mation associated with gas between the nuclei, perhaps in the 
tidal bridge or clump discussed in Section [8] 




Fig. 12. Southern Nucleus: Resulting apparent dispersion of the 
superposition of a 'disc' population and a 'cloud-crushing' pop- 
ulation, both with velocity dispersions of 220 km s" 1 . Left to 
right: disc velocity field, cloud population velocity field, result- 
ing apparent dispersion. Black circles indicate BH positions. The 
right panel shows that a simple 'disc plus cloud' model can pro- 
duce an increase in dispersion, characteristic of that observed 
(Fig. [4]). This indicates that the high dispersion in the northern 
half of the southern nucleus may be an artifact due to the pres- 
ence of two kinematically different stellar populations with sim- 
ilar luminosities. 



7. Jeans Modelling, Nuclear Masses, and Mass to 
Light Ratios 

Knowledge of the rotational velocity field allows one to make 
simple Keplerian mass estimates. However, here the compara- 
tively large velocity dispersions (generally cr/v rot > 1) necessi- 
tate a more sophisticated approach so as not to underestimate the 
actual mass. Two main avenues can be taken: Schwarzschild or- 
bit superposition and Jeans modelling. The former proceeds by 
reproducing the velocity and dispersion maps from a large suite 
of stellar orbits calculated for a given potential, which is varied 
until a match to the observations is achieved. The latter employs 
the Jeans equations, derived by taking moments of the collision- 
less Boltzmann equation. For our data, Jeans modelling is better 
suited to recover the mass distribution, as the accuracy and res- 
olution of our kinematic data is insufficient for Schwarzschild 
modelling. 

A good overview of Jeans models can be found e.g. in 
Bin ney & Tremainel {2008). An important aspect of our ap- 
proach is that we separate the gravitational potential O from 
the distribution of the stars, i.e. we allow for non-stellar mass. 
We assume isotropic velocity dispersion, which can be justi- 
fied because the data do not allow us to say anything about the 
form of anisotropy, and this choice represe nts the least infor - 
mative, and hence least constraining, option (iDejongheU 19861) . 
Another issue is the shape of the stellar and total mass distri- 
bution. A discy system would be characterised by cr/v rot <?c 1, 
and an oblate spheroidal system, in which the dispersion pro- 
vides non-rotational pressure support thickening the disc, by a 
o~/Vrot £ 1. Fig. @] shows that for both nuclei in NGC 6240 
o~/v rot ~ 1, suggesting that the latter is the more physically re- 
alistic c hoice. However, this is non-trivial to implement (see for 
example Ivan der Marel & van Dokkumll2007l) . Instead we make 
the simplifying assumption that both the stellar distribution and 
potential are spherically symmetric, which will allow us to cap- 
ture the characteristics of the system. It also has the advantage 
of yielding an analytical solution. 

In order to assess the uncertainties of our result due to the 
choice of mass distribution, we also ran a set of models assuming 
an axisymmetric thin disc with an isotropic pressure component; 
the enclosed masses found with this model, which can be con- 



10 



H. Engel et al.: NGC 6240: merger-induced star formation & gas dynamics 



sidered to be at the other end of the range of physically plausi- 
ble models, agreed with the results of th e spherically symmetric 
mode ls within a factor of less than two. iHaring-Neumayer et al.l 
(2006) model the gas kinematics of Cen A with three differ- 
ent Jeans models: an axisymmetric thin disc with and without 
pressure terms, and a spherically symmetric model. They find 
that the inferred black hole masses span less than an order of 
magnitude. Comparing their results with those of others, these 
authors find that, whilst favouring the pressure-supported thin 
disc model as the most physically plausible, the spherical Jeans 
model agrees best with t he results from Sch warzschild orbit su- 
perposition modelling by lSilge et al.l (120051) . The uncertainty of 
Jeans modelling results due to the unknown velocity disper- 
sion anisotropy is more difficult to assess due to the degener- 
acy between the int egrate d mass and the anisotr opy p arameter 
p. Both lWolf etaD (l2010h and|Mamon & Boud (l2010h attempt 
to quantify this (see Fig. 1 in IWolf et a0 1201 Ol and Fig. 1 in 
Mamon & Boue 2010) in the case of mass measurements for the 
Carina dSph and DM haloes, respectively, finding that the range 
spanned by results for the integrated mass assuming extreme val- 
ues of j3 are large at small physical radii, decreasing to less than 
a factor of two at radii comparable to or larger than the half-light 
radius. Since a spheroidal system like the one we modelled is 
unlikely to have extremely anisotropic velocity dispersions, and 
since also we are measuring mass at radii larger than the half- 
light radius, we can assume the factor of two derived by these 
authors to be an upper limit on the error due to our choice of 
f3=0. We thus estimate that the uncertainties introduced by the 
assumptions inherent in our modelling are unlikely to be larger 
than a factor of two. 

Jeans modelling requires luminosity, rotational velocity, and 
velocity dispersion profiles as input. Since the region between 
the two nuclei is likely strongly perturbed, we measured az- 
imuthally averaged profiles using only a 1 80° wedge in the out- 
side halves (i.e. opposite the merger centre) of the galaxy, after 
subtracting the continuum associated with the recent starburst 
(assuming an intrinsic Wbtj of 22 A, as measured at the knot of 
Bry emission north-west of the northern nucleus, cf . ^9]) and cor- 
recting for extinction. To these we fitted a Sersic function to the 
profiles, constraining it with the NIRC2 profile at > 0.4", and 
with the extinction-corrected SINFONI profile at < 1.0" (i.e. 
some overlap between the two profiles is included). As part of 
the fitting process, the model profile was convolved in 2D with 
the SINFONI PSF, to account for beam smearing at small scales. 
We then analytically deprojected these LOS luminosity profiles. 
The measured rotation curves and dispersion profiles need to be 
corrected for beam smearing and projection along the line of 
sight in order to recover the intrinsic kinem atics. For this, we 
used the code described in lCresci et al.l (l2009h . 

With these inputs, we compute M(r). The total mass en- 
closed within the cut-off radius (250 pc and 320 pc for northern 
and southern nuclei, respectively) is found to be 2.5 x 10 9 M© and 
1.3 x 10 10 M o . This is comparable to that se en in the central few 
hundred parsecs of nearby AGN (Figure 7,|Davieset ID[2007). 
We also calculate a global K-band mass-to-light ratio by find- 
ing the multiplication factor that best matches the luminosity en- 
closed at radius r to M(r). Here, as in the rest of the paper, L K is 
taken to be the total luminosity in the 1 .9-2.5/mi band in units of 
bolometric solar luminosity where 1 L = 3.8xl0 26 W (we note 
that a commonly used alternative is the definition via solar K- 
band lumnosity density, 2.15xl0 25 W jim~ l ). We obtain values of 
5.0 Mq/Lq and 1.9 M /L Q for the northern and southern nuclei. 
iDownes & Solomonl {1998) find typical gas fractions of -15% 
for local ULIRGs; in conjunction with these modelling results, 



this implies total stellar masses of 2. 1 x 10 9 M o and 1 . 1 x 10 10 M o , 
and stellar mass-to-light ratios of 4.3 M /L o and 1.6 M o /L for 
the northern and southern nuclei, respectively. We discuss the 
possible inferences from these results regarding the nature of the 
nuclei and the progenitors' Hubble types in CT01 

It is tempting to compare the mass enclosed in the inner- 
most few parsecs to the expected black hole masses. One av- 
enue to estim ate the black h ole m asses is afforded by the X-ray 
luminosities. Vig nati etaD (Il999l) measure a total absorption- 
corrected nuclear X-ray luminosity in the 2-10keV ran ge of 
3.6xl0 44 ergs _1 with BeppoSAX. iKomossa etaD (l2003l) . us- 
ing CHANDRA, find that the northern and southern nuclei 
have absorption-corrected 0.1-10keV X-ray luminosities of 
0.7xl0 42 erg s" 1 and 1.9xl0 42 erg s" 1 , respectively. The discrep- 
ancy is most likely due to different absorptio n corrections 
(IKomossa et al] 120031 derive Nh ~ 10 22 cm~ 2 , I Vignati etaL 
1 19991 measure N H ~ 2 x 10 24 cm" 2 ). Since IVignati et aL 
(1999) derive their absorption correction from a significantly 
larger wavelength range (up to lOOkeV) we use their values 
for the following estimates. We convert the X-ray luminos- 
ity to monochromatic 5 10 A luminosi t y using the luminos- 
ity d ependent a ox -relation (ISteffen et all [20061; iMaiolino et all 
2007), and then convert to AGN bolometric luminosity adopt- 
ing L hol =lvL v (510 0A) (e.g. iNetzer & Trakhtenbroti l2007h . 
Das yra et"aD (l2006bl) calculate Eddington ratios for a sample of 
34 ULIRGs, finding a median value of 0.37. With this value for 
Lboi/LEdd, we arrive at a combined black hole mass of 4xlO 8 M . 
However, the uncertainty on this estimate is at least a factor of a 
few. 

Alternatively estimates can be derived from t he stellar veloc- 
ity d ispersions. Using the M##-cr relation from lTremaine et aD 
( 2002) and the stellar dispersions measured at the locations of 
the AGN (-200 km s" 1 and -220 km s" 1 for the northern and 
southern nuclei, respectively) yields expected central black hole 
masses of 1.4±0.4xl0 8 M o and 2.0±0.4xl0 8 M o for the north- 
ern and southern nuclei, respectively. The uncertainties on black 
hole masses derived from the M##-cr relation are also a factor of 
a few. 

Finally, a measurement (with a rather smaller uncertainty) 
of Mbh m the southern nucleus has recently been made using 
high resolution stellar kinematics. Resolving the gravitational 
sph ere of influence of the black hole and modelling the kinemat- 
ics, |Midiin^itaLl d20lo|) derive a mass of 2.0 ± 0.6 x 10 9 M . 

Although the range of Mbh above spans an order of magni- 
tude, the very large uncertainties of the first 2 estimates means 
that all the estimates are formally consistent (within ~2<x). We 
therefore caution against over-interpretation. 

8. Merger Geometry & Stage, And Cold Gas 
Concentration 

Starting from the assumption that the project ed angular m o menta 
of the nuclei are the true angular momenta, iTecza et 
proposed a merger geometry in which one of the nuclei is copla- 
nar and prograde with the merger orbital plane, and the other is 
inclined with respect to the orbital plane (their Fig. 11). They 
further supported the notion of at least one of the galaxies be- 
ing subject to a prograde encounter by noting that the formation 
of tidal tails such as those of NGC 6240 is favoured in prograde 
encounters. Here we would like to build onto and expand this 
discussion. 

NGC 6240 indeed has long, well de veloped tidal t a ils, as 
can be seen in optical and HI data dGerssen et all 120041; 
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Yun & Hibbard, 2001), covering a projected size of ~45kpc 
from north to south, about half that of the Antennae. In the HST 
images, a continuous dust distribution across the centre strongly 
suggests that what we are seeing is an extended, well devel- 
oped single tail curving in front of the system, rather than a 
series of disconnected shorter features. However, NGC 6240 's 
tails are not as dense or long as those of e.g. the Antennae. 
What can we derive from this with regard to the merger ge- 
ometry? It is generally known that in merger simulations, the 
more coplanar/prograde, and thus 'resonant', an encounter is, 
the stronger the resulting tidal tails are. It is quite clear however, 
from looking at the relative orientation of the rotation axes of 
the nuclei, that NGC 6240 cannot be a perfectly coplanar merger. 
But a merger must not be exactly coplanar and prograde in or- 
der to produce tidal tails - an example of this is NGC 7252, 
which has tidal tails akin to those of NGC 6240. NGC 7252's 
tidal featur es were successfully reprod uced through N-body sim- 
ulations by Hibbard & MihosJ (1 19951) ; their simulation had one 
of the progenitors on a coplanar, prograde encounter, but the 
other galaxy at an inclination of ~ 45 degrees to the orbital plane 
- this illustrates that there is a non-negligible merger parame- 
ter space around 'coplanar/prograde' which will also result in 
pronounced tidal tails. Therefore, whilst not pinning down the 
merger geometry exactly, the existence of such long and well 
developed tidal tails is convincing evidence that the merger ge- 
ometry must be within a reasonable parameter space around pro- 
grade/coplanar. The stellar kinematics at first sight seem to dis- 
agree with this, giving the appearance that the nuclei are counter- 
rotating. However, co-rotating nuclei would also look like this if 
one nucleus is inclined behind the plane of the sky, and the other 
in front. Thus the stellar kinematics are consistent with the view 
that NGC 6240 is not too far from being a prograde merger. That 
the nuclei are clearly not exactly aligned could be the reason 
for the additional shorter tails that make the system look rather 
messy. 

We therefore conclude that NGC 6240 's merger geometry 
most likely tends towards coplanar/prograde, and in this sense 
is perhaps quite simil ar to the geometry proposed for NGC 7252 
bv lHibbard & Mihosl (119951) . 

It is more difficult to deduce the stage of merging the two 
galaxies comprising NGC 6240 are currently in - we cannot pin 
down their merger geometry exactly, as discussed above, and 
neither do we know the 3D relative velocities of the two nuclei, 
so we cannot calculate whether the two nuclei will merge imme- 
diately or separate again before eventually coalescing. We know 
that they must certainly have already experienced their first close 
encounter, since tidal tails only form after the first strong gravi- 
tational interaction of the two galaxies. This is also supported by 
the high luminosity of the system, because elevated levels of star 
formation also are only expected after the first encounter leads 
to collision and compression of gas complexes. And the fact that 
the two nuclei are still separated by about a kpc in projection in- 
dicates that the system is not yet at the 'final coalescence' phase. 

The observe d peak of the gas emission between the nuclei 
(Fig. [5]), which pLicconi et all (Il999l) found to display a velocity 
gradient and interpreted as due to a self-gravitating gas disc lo- 
cated in between the nuclei, is puzzling, since from merger sim- 
ulations this is generally not expected to occur. We could be see- 
ing a collapsed gas clump in a tidal arm - such features are found 
to occur quite often in simulations of merging galaxies above 
a cer tain gas fraction dWetzstein et all 2007; Bournaud et all 
2008), and are also seen in observations (Knierman et al., 2003). 
An alternative, and perhaps more likely, interpretation of this gas 
concentration is a tidal bridge connecting the two nuclei that is 



viewed in projection. Such tidal bridges are often seen in merger 
simulations. However, simulations have never yielded a system 
in which the total gas mass is dominated by that between the nu- 
clei. And although bridges are produced most strongly in a per- 
fectly prograde coplanar encounter (which, as discussed above, 
is inconsistent with our observations), as one deviates from a 
prograde coplanar geometry the strength of the bridge decreases. 

An alternative, although somewhat rarer, way to drive a sig- 
nificant gas mass away from the nuclei is a direct interaction 
(i.e. low impact parameter), in which the ISM of the progenitors 
collides. The best example o f this is UGC 12914/5 (the Taffy 
Galaxies; iBraine et all l2QQ3h . in which about 25% of the CO 
emission originates between the nucle i. In this system, t he pro- 
genitors collided at about 600 km s" 1 . IBraine et al l(l2003l) argued 
that while this would have ionised the gas, the cooling time is 
short enough that H2 could reform while the discs are still pass- 
ing through each other. From the ratio of 12 CO to 13 CO these 
authors found that the optical depth of CO in the bridge was far 
less than on the nuclei. Combined with the effect of increased 
CO abundance due to grain destruction, this led them to suggest 
that the actual mass of molecular gas may be rather less than im- 
plied by the CO luminosity, perhaps only ~6% of the total gas 
mass. Since direct collisions such as these tend to create ring 
galaxies (or in the case of the Taffy Galaxies, an incomplete ring 
or hook), this scenario seems unlikely for NGC 6240. 

In NGC 6240, the CO map indicates that CO luminosity is 
dominated by that between the nuclei. We must therefore con- 
sider physical effects that could significantly affect the flux-to- 
mass conversion factor between the two regions. It is already 
known that sub-thermal emission from non-virialised clouds can 
radically m odify the CO- t o-H2 conversion factor, as shown in 
Fig. 10 of iTacconi et al.l (120081) . Such an effect could easily 
boost the emission from gas in a drawn-out, thin and diffuse gas 
bridge compared to denser material around the nuclei. And it is 
supported by other observations: a comparison of the CO line 
emission (Fig [5]) w ith the 1.315 mm continuum measured by 
ITacconi et alld 1999b (their Fig. 3, left panel) shows that the dust, 
as traced by the mm continuum, is concentrated on the nuclei 
rather than following the CO line emissi on. Furthermore, in a de- 
tailed analysis of various CO transitions. iGreve etal (2009) were 
unable to find a single set of average H2 conditions that comes 
close to reproducing the observed line ratios. A second impor- 
tant issue is the impact of X-ray irradiation on the CO abundance 
arou nd the nuclei, since the AGN are quite luminous. iLutz et aD 
(12003b argue that as little as 25-50% of L bo i is due to the AGN. If 
we take the lower end of this range, and distribute it equally be- 
tween the two AGN, we find each radiates a t L ~ 3 x 10 44 erg s" 1 . 
We also adopt a column of 2 x 10 24 cm -2 dVignati et all Il999b 
and assume a gas density of 10 5 cm -3 . Under these conditions, 
hard X-rays from the AGN will cause sufficient ionisation to re- 
duce the CO abundance an order of magnitud e below its typical 
value of 10" 4 out to a radius of 250 pc (see lBoger & Sternbergl 
I2005L and Davies & Sternberg in prep.). Thus, even if there is 
significant gas mass around the nuclei (as the models imply), 
one would expect rather little CO emission. That these effects 
will conspire to make the CO luminosity distribution look rather 
different to the molecular gas mass distribution should be borne 
in mind. 

9. Scale of Starburst & Star Formation History 

While lTecza etaD (12000b concluded that the nuclei were the pro- 
genitor bulges, they did not quantify how much of the K-band lu- 
minosity arises in the recent starburst and how much is due to the 
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progenitor bulges themselves. Our data enable us to resolve this 
issue. Our method is to measure an average value for Wsry away 
from the nucleus where dilution from the bulge is small, and 
assert that this is representative of the whole starburst. The fun- 
damental assumption is that, integrated over a sufficiently large 
aperture, we are summing a fair cross-section of star clusters 
and hence probing the average star formation properties. This 
minimises the impact of the stochastic nature of individual clus- 
ters. We have therefore chosen a region containing a high density 
of clusters, and measured W# rr in a 0.8" aperture correspond- 
ing to 400 pc so that many clusters are included. Such a region 
lies a bout 1" to the west of the northern nucleus (Poll ack et all 
120071) . sufficiently far that the K-band continuum from the nu- 
cleus is very faint. Fig. [13] shows that, despite the Bry flux being 
higher on the nuclei, both here (and also in other off-nuclear re- 
gions) Wbtj = 22 ± 6 A. We use it to obtain the fractional contri- 
bution of the young stellar population in 1" diameter apertures 
centred on the nuclei, by dividing the measured W# rr by this in- 
trinsic value. We then find the K-band luminosity emitted by the 
starburst population to be 3.6 x 1O 8 L (northern nucleus) and 
1.3 x 1O 9 L (southern nucleus), implying L K y 0ung /L K j ota i to be 
0.36 (northern nucleus) and 0.32 (southern nucleus) - i.e. only 
about 1/3 of the total K-band luminosity of the nuclei is due to 
recent star formation. 

This implies that the K-band luminosity of the nuclei is dom- 
inated by a population of stars older than 20 Myr. An important 
question is whether this is consistent with claims that one needs 
supergiant templates in order to match the CO bandhea d depth 
(Lester et allll988USugai e7aflll997l:lTecza et alll2000b . In the 
following, we show that basing conclusions about the dominant 
stellar popula tion on a fit to a sin gle t emplate is very uncertain . 
The tables in lOriglia et all (1 19931) and lForster Schreiberl (l200Qb 
show that the equivalent widths of the CO bandheads vary con- 
siderably between individual stars; and it is also well known 
that higher metallicity populations have deeper bandheads be- 
cause the K-band c ontinuum is domina ted by cooler sta r s. The 
templates found by iTecza et al.1 (l2000b and ISugai et al.l (1 19971) 
to provide the best fit to the K-band absorption features are of 
K4.5Ib and K2.5Ib stars, respectively: i.e. specifically K type 
sub-luminous supergiants. Both authors infer from this that the 
domin ant stellar population is late type supergiants. ITecza et al.l 
(2000) go on to conclude that these were produced as a re- 
sult of a burst of star formation that lasted for 5 Myr and oc- 
curred ~20Myr ago. We have examined the Tecz a et al ] (120001) 
starburst scenario in detail using the population synthesis code 
STARS, and found that it would lead to ~ 60% of the K-band 
luminosity being due to M supergiants (cooler than 4000 K) and 
only ~3% due to K supergiants (temperatures 4000-4600 K). 
Thus the scenario of a ~ 20 Myr old starburst leads to two ap- 
parent contradictions: How can M supergiants dominate the lu- 
minosity and yet have absorption features too deep to match the 
spectrum? And how can K supergiants provide a good match to 
the spectrum and yet contribute only an insignificant fraction of 
the luminosity? The resolution is simply that the galaxy spec- 
trum consists of many different types of stars, and a single tem- 
plate can at best be characteristic of the sum of all these. The 
important conclusion here is that while the composite spectrum 
is best matched by a supergiant template, other types of stars still 
make up a significant fraction of the near-infrared continuum. 

In addition to the current scale of the starburst, another key 
point of interes t are the pas t an d future st a r formation rates. 
iDi Matteo etaD (120071 120081) and lCox etal] (I2006L 120081) offer 
a large sample of medium-resolution simulations comprehen- 
sively covering a wide range of gas fractions, merger geometries, 



and mass ratios, as well as employing different codes. Their syn- 
thesised results regarding the qualitative evolution of the star for- 
mation during a merger are that nearly all encounters are marked 
by at least two peaks in the star formation rate; one at the first 
encounter, and one upon final coalescence. However, the relative 
strength of the two peaks differs, depending on th e merger geom- 
etry; but on average the second peak is stronger ( Di Ma tteo et"ail 
2007, 2008). This is particularly evident for mergers that tend 
towards coplanar/prograde geometries, as appears to be the case 
for NGC6240. 

The aforementioned simulations were carried out at medium 
resolution, and thus one might ask whether their results would 
change when going to higher reso lution simulations res olving 
the multiphas e ISM . However, as Bou rnaud et al.l (|2008) and 
iTeyssier et al.l (l2010h show, the qualitative evolution of the star 
formation rate during a merger remains largely unchanged, the 
only significant difference is that star formation proceeds much 
more effectively, resulting in predicted star formation rates a 
factor of up to ten times larger than those predicted by lower- 
resolution simulations (iTeyssier et all l2010b . 

We therefore conclude that NGC 6240, being between first 
encounter and final coalescence, most likely experienced its first 
peak in star formation rate in the recent past (triggered by the 
first encounter); has currently elevated levels of star formation 
compared to a quiescent galaxy; and will experience another, 
likely stronger, peak in star formation rate in the near future 
when the galaxies coalesce. This is supported by the observed 
Bry emission, which indicates that star formation must cur- 
rently still be on-go ing, and by measurements of cluster ages 
(iPollack et all 12007) which are found to be typically very young 
- in a population of clusters with a range of ages, it is the brighter 
ones that are more easily detected, and because clusters fade 
quickly, these will also inevitably be younger. We furthermore 
note that since NGC 6240 is already just below the canonical 
ULIRG threshold of L IR > 10 12 L o , it is safe to predict that it will 
breach this threshold once the final starburst is triggered, and 
will become a bona fide ULIRG. 

In order to be able to make more quantitative analyses, 
we use STARS to calculate a range of observables from two 
star formation histories, which both display the generic fea- 
tures discussed above; an initial peak at first encounter fol- 
lowed by a gradual rise or plateau, and a final intense burst. 
One of these simulations is selected from a set of simulations 
intended to reproduce t he properties of the Antennae galaxies 
dKarl et all 120081 l2010h . and the other one from the library of 
Joh ansson et al.l (120091) . Neither are coplanar/prograde, although 
they tend towards it, and both produce well-developed tidal tails 
during the interaction. The rationale of using star formation his- 
tories from two different simulations is to provide an estimate of 
the uncertainties introduced by the quantitative differences in- 
troduced by the exact simulation details. Since we only know 
that NGC 6240 must be between first encounter and final coa- 
lescence, we use the full range of simulated properties between 
these two points, as well as our measurement and associated un- 
certainties of W# rr , as constraints to derive the uncertainties in 
the star formation properties. 

In Table [2 we have used these star formation histories to cal- 
culate the mass of young stars formed in both scenarios, as well 
as the current star formation and supernova rates. To do so, we 
have applied a scaling so that the K-band luminosity matches 
that observed for each nucleus, under the assumption that the 
evolution of the central SFR mirrors that of the global SFR. The 
table shows that quantitatively the results for the two different 
star formation histories do not differ greatly - in fact, the varia- 
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tion within each scenario over the possible time range is larger 
than the variation between the averages of the two cases. 

One conclusion is that the bolometric luminosity of the star- 
burst from the two nuclei together is only 20-30% o f the sys- 
tem's total. This follows in the same direction as iLutz et all 
( 2003) that the recent starburst only contributes part of the to- 
tal luminosity. W e estimate a somewhat lower fraction than 
ILutz et all (120031) . who constrain starburst contribution to be 50 
to 70%, possibly due to the fact that we are looking at the central 
region whereas they investigated the global luminosities. 

An important co ncern is that our result s appear to be incon- 
sistent with those of iBeswick et al.l ([2001) who find, based on 
the 1 .4 GHz continuum, significantly larger star formation and 
supernova rates of 83.1 M© yr" 1 and 1.33 yr" 1 . We argue that 
this is due to the different star formation histories adopted. We 
show that this can have a major impact on interpretation of the 
data, and that accounting for it makes the radio continuum data 
consistent with our results. 

We first address a minor correction, specifically a ~10% 
increase in radio fl ux due to the larger aperture used by 
IBeswick et al.l (1200 ll) . which included a substantial flux contri- 
bution from an off-nuclear source not included in our apertures 
('N2' in their nomenclature, see their Fig. 2). A lar ger correction 
may b e needed by analogy to Arp220 for which iRovilos et all 
( 2005) examined the 18-cm lightcurves of the supernovae. They 
found the type Iln RSNe model insufficient and arrived at a su- 
pernova rate that was a factor of ~3 smaller than the ~2yr _1 
often quoted. A similar correction may be applicable to the rates 
inferred from radio measurements for NGC 6240. This may be 
related to the third issue, wh ich concerns the star formation 
history. The formulae used by IBeswick et al.l (1200 ll) to convert 
1 .4 GHz lum inosity to star fo rmatio n and sup ernova rates (based 
on work by ICondon & Y in 1990; ICondon|[l992l: ICramlll998l) 
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were derived empirically based on measurements of normal disc 
galaxies for which a constant star formation rate is characteristic. 
This implicit constant star formation rate is very different to the 
increasing star formation rate found in merger simulations. We 
show below that this can make as much as a factor of three dif- 
ference in the star formation and supernova rates that are derived 
from observables. 

We use STARS to estimate the star formation rate from a 
given K-band luminosit y for three dif ferent star formation his- 
tories. First, following IBeswick et al.l (120011) . we adopt contin- 
uous star formation for a durat i on of 20 Myr, correspondin g 
to that estima ted by iTecza et al . (2000), Pasquali efalj (|2003|), 
iPollack et al.l (120071) for the most recent star burst. In order to 
reach L^=1.7xlO 9 L , corresponding to the total measured for 
the starburst population, we require a constant SFR of 25 M 
yr -1 . On the other hand, for the recently increasing star forma- 
tion rate typical of the merger scenarios, we find star formation 
rates of ~ 10M o yr -1 . Thus we find a factor ~ 2.5 difference in 
the SFR required to reach the same K-band luminosity, depend- 
ing on the star formation history. We can estimate the supernova 
rates in a similar way: For 20 Myr of continuous star formation at 
25 M yr -1 , we calculate a supernova rate of 0.3 yr -1 . In contrast, 
for our merger scenarios, we find supernova rates of ~ 0. 13 yr -1 . 
Thus, much of the discrepan cy between our supe rnova rate and 
the (corrected) value from IBeswick et al.l (120011) is due to the 
different star formation histories adopted. 

This illustrates an important caveat: the choice of star forma- 
tion history can have a far-reaching impact on the interpretation 
of the observables. We note that the very similar derived stellar 
population parameters for our two different merger star forma- 
tion histories indicate that, for mergers, using a 'characteristic' 
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Fig. 13. Left: Map of Bry flux (units given in 10 16 W m 2 mum). 
Right: Map of W# rr (units given in A). White contours trac- 
ing the K-band stellar continuum, and black circles denoting the 
black hole locations are overdrawn on both maps. Although the 
most intense Bry emission is on the nuclei, the equivalent width 
here is lowest, suggesting that the K-band light is dominated by 
an older stellar population. We note that in §0 we showed that 
dilution by hot dust emission associated with the AGN cannot 
play a role in reducing Wg rr . 
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Fig. 14. Spectrum of the region to the northwest where the Wsry 
is not diluted by continuum emission from the progenitors, in- 
tegrated over an 0.8" diameter aperture. The flux is in units of 
10" 16 W m~ 2 mum~ l . The measured Wsry is 22 ± 6 A. 



star formation history (i.e. an initial peak at first encounter fol- 
lowed by a gradual rise or plateau, and a final intense burst) will 
yield much more reliable quantities than either instantaneous or 
continuous star formation models. 



10. Stellar Populations of NGC 6240 

Did the Progenitors Have Bulges? The existence or otherwise 
of bulges (where 'bulge' here includes both pseudobulges and 
classical bulges) in the galaxies und ergoing a merger ha s a pro - 
found impact on the star formation (Mi hos & Hernquisl [T996). 
Spri ngel et al.l (I2005I) show in their Fig. 15 that in a major merger 
of discs (gas fraction 10%) without bulges, the first encounter 
already triggers an episode of enhanced star formation that lasts 
for 300-400 Myr. The second encounter triggers a more intense, 
but also somewhat shorter (-50 Myr), starburst 800-900 Myr 
later. In a simulation that included a bulge and a central black 
hole, there was almost no perceptible change in the star forma- 
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Table 2. Characteristics of Starburst Populations 



simulation nucleus L Kyoung L bo i young mass of young stars SFR SNR 

[Lp] [L q ] [M ] [Mayr- 1 ] [yr' 1 ] 

"^T North 3.6x10* 2.9^ x 10 1U x 10* 23™ 0.03+0.1 

'2' North 3.6xl0 8 3.1+^ x 10 10 3.7+0.3 xlO 8 2.2+°* 0.03+0.01 

'1' South 1.3xl0 9 10.0+2.0 xlO 10 1.5 + g-f x 10 9 8.2 + ^ 0.11+°,$ 

'2' South 1.3X10 9 ll.O+^xlO 10 1.3+0.1 XlO 9 7.8+j;? O-lO^m 

References. L Kjoung was measured in 1" diameter apertures using the dereddened continuum flux map and adopting an intrinsic W Bry = 22 A (see 
|9]). Values for L boh stellar mass, and SFR were inferred from this together with the star formation history from the simulations at all points 
between first encounter and coalescence; given are the median value and the range sampled. 

Table 3. Properties of the Nuclei 



nucleus 


^max 


CT 


n 


6 


r eff 


Thin 


L K 


M dyn 


M dyn /L K 




[km s" 1 ] 


[kms" 1 ] 






[pc] 


[pc] 


[L Q ] 


[L Q ] 


[Lq/Mq] 


North 


205 


200 


2.8 


0.3 


490 


265 


l.lxl0 y 


2.5xlO y 


5.0 


South 


300 


220 


2.8 


0.35 


180 


290 


5.9xl0 9 


1.3xl0 10 


1.9 



References. v max is the peak measured stellar velocity, cr is the average/typical stellar dispersion ( 32.1.41 ); n and r ef f are the Sersic index and 
effective radius of the light profile (§0, e is the ellipticity (|5j), r kin is the extent to which rotation is seen, L K the total K-band luminosity (©, and 
Mdyn the dynamical mass as derived from the Jeans modelling (§|7]). 



tion rate at the first encounter. This remained constant at ~5- 
10% of the peak star formation rate until the second encounter, 
when the starburst is triggered and the peak star formation rate 
was achieved. Because there had been no earlier phase of en- 
hanced star formation, this starburst is more intense than it would 
be without a bulge. Which of these two different scenarios holds 
true for NGC 6240? 

The evidence overwhelmingly points towards the progeni- 
tors having had a prominent bulge component. First, the majority 
of massive spirals at low redshifts have bulges, although at z > 2 
the situation may well be different as bulges may not yet have 
had time to grow (iGenzel et all 120081) . Second, we know that 
both nuclei host an AGN and therefore harbour a supermassive 
black hole. The implication of the Mbh~o~* relation is therefore 
that the progenitors must have had bulges. Using estimates of the 
black hole masses based on the AGN luminosity and our mea- 
sured stellar dispersion (©in conjunctio n with the empirica l 
correlation between M B h and bulge mass ( Haring "&Ril 120041) 
yields associated bulge masses of 8.7xl0 10 M o and 1.2xlO u M 
for the northern and southern nuclei respectively, an order of 
magnitude larger than their dynamical masses. However, due to 
their transformational impact on both black hole and host galaxy, 
it is uncertain whether such scaling relations can be expected to 
hold during on-going mergers. An argument favouring the exis- 
tence of bulges in the progenitors stems from the merger simu- 
lation: the galaxies in our simulation had bulges, and the good 
match to the observe d Ws r y suggests that th e st ar formation his- 
tory is reason able. As lSpringel et al.l (120051) and lDi Matteo et all 
(120071 120081) show, bulgeless mergers have markedly differ- 
ent star formation histories, with the first encounter produc- 
ing the most potent burst of star formation (it should further- 
more be noted that in higher-resolution simulations, the star- 
burs t triggered by the fi nal coalescence occurs somewhat earlier; 
cf. iTeyssier et al.ll2010h . Thus, at the current phase, we would 
expect to see significantly lower W# rr if the first encounter had 
indeed triggered the more powerful star formation event. 

Are the Nuclei the Progenitors' Bulges? Are the nuclei we see 
(remnants of) the progenitor galaxies' bulges, did they originate 



only in the first encounter a few lOOMyr ago, or is the rota- 
tion we see the signature of stellar discs? The latter is unlikely, 
since the high dispersion clearly indicates a thickened structure, 
and nuclear stellar discs in contr a st are often accom panied by 
sigma-drops (Pel etier et alll2007l) . lTecza et al.l (l200Qh proposed 
the first, based on an estimate of the global mass-to-light ratio of 
the nuclei. Our high spatial resolution data allow us to explore 
this in more detail using not only the mass and luminosity, but 
also the effective radius and kinematics, and our knowledge of 
the stellar populations. 

Table [3] quantifies the kinematics, luminosities, masses, and 
sizes of the two nuclei. How do these compare with the charac- 
teristics expected of a bulge? An argument in favour of the nuclei 
being bulges is supplied by estimates of the mass-to-light ratio 
of the non- starburst stellar population: correcting the global stel- 
lar mass-to-light ratios (§|7] assuming a gas fraction of -15%) 
for the starburst contribution (Table |2]), we find that the non- 
starburst nuclear populations must have mass-to-light ratios of 
~ 5.2 Mq/Lq (northern nucleus) and 1.8M /L© (southern nu- 
cleus) - noting that the associated uncertainties are estimated 
to be about a factor of two. If the nuclei only formed after the 
first encounter (~560Myr ago in the 'Antennae' scenario), they 
should have a M/L K of at most 2M O /L , whereas a several Gyr 
ol d stellar populatio n has a M/L K of ~6-lOM /L o (e.g. Fig 4 
in iDavies et al.l 120071) . This, as well as the small relative con- 
tributions of the starburst population to the nuclear masses and 
luminosities, indicates that a substantial fraction of the nuclei 
consists of stars predating the beginning of the merger. 

NGC 6240 has often been cited as having anomalously 
high stellar dispersion. However, excluding the region towards 
the northern half of the southern nucleus which we discuss 
in §0 the dispersions typical for the nuclei are -200 km s" 1 . 
These are not unusua l for late-type bulges, as Fig. 18 in 
Kor mendy & Kennicut t ( 2004) shows. Placing the two nuc lei on 
the 6 vs. v max /cr diagram of lKormendy & Kennicuttl (120041) (their 
Fig. 17), it appears that the northern nucleus lies comfortably 
within the region typical for late-type pseudobulges, whereas 
the southern nucleus would represent a somewhat more extreme 
case, having a higher v max /cr than other galaxies. In both cases, 
however, the kinematics strongly support the case for a disc- 
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like, secularly-grown pseudobulge rather than a merger-formed 
classical bulge. The Sersic indices are rather typical of classi- 
cal bulges, but we caution that the luminosity profiles are dis- 
turbed and influenced by recent star formation, and hence the 
values derived from the Sersic fits should be treated with cau- 
tion. For the same reason, we prefer to measure the sizes from 
the extents to which rotation can be seen (indicating the extent 
to which the nuclei dominate the light distribution), rather than 
the effective radii obtained from the Sersic fits; this gives radii 
of 265 pc and 290 pc for the northern and southern nuclei, re- 
spectively (as opposed to 490 pc and 1 80 pc based on the Sersic 
fits). Thi s places them at the lo wer end but within the range 
given by iMoriondo et al.l (Il999l) who measure the bulge radii 
for 40 spiral galaxies and find values ranging from 230 pc to 
1.93 kpc with an average of 650 pc. The same authors also look 
at the bulge masses, where they find values from 1.4xl0 9 M o to 

I. 8xl0 n M o with an average of 2.5xlO lo M . Our analysis ($7^ 
yields masses of 2.5xl0 9 M o and 1 9xl0 1Q M o , aga in at the lower 
end but within the range found by Moriondoetal] Jl999). 

Would we expect the bulges to survive to this merger state? 
Comparatively little is known about t he evolution of galaxy 
bulges during mergers. Observationally. lAguero et all (1200 ll) re- 
port velocity dispersions of 217 km s" 1 and 280 k m s" 1 in the 
two nu clei of the merger system AM 2049-691, and I Shier et al.l 
(1998) investigate 11 multiple-nucleus merging objects and find 
nucleus sizes between 200 and 1600pc (on average 600 pc), 
velocity dispersions of 66-151 km s" 1 (excluding NGC 6240, 
which was also sampled) with an average of 109 km s" 1 . 
However , these w e re to a large part galaxies in earlier merger 
stages. Hernquist (1993) simulates the merger of two spiral 
galaxies with bulges. His model suggests that the bulges remain 
very nearly intact until final coalescence. The spin (rotational 
velocity) is reported to remain nearly constant throughout the 
simulated merger. 

In light of this and the fact that the nuclei are found to- 
wards the lower end of the range of masses and sizes found for 
a large sample of spiral galaxy bulges, we conjecture that the 
nuclei are indeed the progenitor galaxies' bulges, most likely 
pseudo-bulges, that have had their outer, less tightly bound lay- 
ers stripped off during the merger, leading to the sub-average 
sizes. Dispersion may have increased due to 'heating' during the 
interaction by gravitational impulse, and a gravitational poten- 
tial deepened by influx of gas into the central region. Since the 
progenitors did not only have a prominent bulge component, but 
must also have had substantial discs in order to allow such pro- 
nounced tidal features as we observe to develop, the progenitors 
were most likely of Hubble type Sa/Sb. 

II. Conclusions 

We present new adaptive optics integral field spectroscopy near- 
IR data and CO(2-l) interferometric line observations of the 
nearby merger system NGC 6240. Our main conclusions are: 

- The ICalzetti et al.l (l200Qh reddening law provides the best 
fit to photometric data points spanning 0.45/mi to 2.22/zm. 
The spatially resolved extinction is generally moderate, with 
A K ~ 2-4. 

- The locations of the stellar kinematic ce ntres are consisten t 
with the black hole locations proposed by Max et al. (2007). 
However, an additional stellar population along the line of 
sight may be perturbing the velocity field of the southern 
nucleus and causing the higher dispersion across the northern 
side. 



- Jeans modelling of the stellar kinematics, assuming spherical 
symmetry and isotropic dispersion, gives M = 2.5 x 10 9 M o 
and M/L K = 5.0 M o /L for the northern nucleus (out to 
250pc), and M = 1.9 X 10 10 M o and M/L K = 5.0 M /L o 
for the southern nucleus (out to 320 pc). 

- The presence of tidal arms, and the still- separated nuclei in- 
dicate that NGC 6240 must be between first encounter and 
final coalescence. The stellar velocity field, and strength and 
prominence of the tidal arms indicate that this system is not 
perfectly, but reasonably close to prograde coplanar. 

- Although simulations can produce significant projected gas 
mass between the nuclei, these are never as prominent as the 
observed concentration in CO luminosity. In order to explain 
this, additional physical effects such as spatial differences in 
CO-to-H2 conversion factor or CO abundance are needed. 

- Star formation histories from numerical simulations display 
generic features. Using these in combination with our con- 
straints on the merger stage and measurements of Wsry, we 
calculate the properties of the young starburst population. 
We find that recent star formation accounts for only about 
1/3 of the total K-band luminosity. Thus the stellar luminos- 
ity is dominated by stars predating the merger. 

- Whereas the differences between the starburst properties in- 
ferred from different merger star formation histories appro- 
priate to NGC 6240 are markedly small, adopting a constant 
star formation rate yields very different results (such as a ac- 
tor 3 change in the derived SFR). We conclude that, when 
characterising the star formation properties of a merger, a 
'generic' merger star formation history should be adopted in 
preference to either instantaneous or constant star formation 
histories. 

- After accounting for the recent star formation, the mass-to- 
light ratios of the remaining stellar population are 5.2 M© /L© 
and l.SM Q /L Q for the northern and southern nuclei respec- 
tively. This implies that a population of stars older than 
~1 Gyr contributes the majority of the nuclear stellar masses 
and luminosities. Combined with the measured size and V/cr, 
as well as results from the simulations, it implies that the two 
nuclei are the remnants of the bulges of the progenitor galax- 
ies. 

Acknowledgements. The authors would like to thank the anonymous referee for 
a thorough and thoughtful reading of the manuscript, which helped to clarify 
and improve the paper. H. Engel would like to thank Scott Tremaine and Pay el 
Das for interesting and helpful discussions on Jeans modelling. The numerical 
simulations were performed on the local SGI-Altix 3700 Bx2, which was partly 
funded by the Cluster of Excellence: "Origin and Structure of the Universe". 
Some of the data presented here were obtained at the W. M. Keck Observatory, 
which is operated as a scientific partnership among the California Institute of 
Technology, the University of California, and NASA. This work was supported 
in part by the NSF Science and Technology Center for Adaptive Optics, managed 
by the University of California at Santa Cruz (under cooperative agreement no. 
AST-9876783). 

References 

Abuter R., Schreiber J., Eisenhauer E, Ott T., Horrobin M., Gillessen S., 2006, 

NewAR, 50, 398 
Agiiero E. L., Diaz R. J., Paolantonio S., 2001, PASP, 113, 1515 
Alonso-Herrero A., et al., 2006, ApJ, 650, 835 
Armus L., Heckman T. M., Miley G. K., 1990, ApJ, 364, 471 
Armus L., et al., 2006, ApJ, 640, 204 
Barnes J., 2002, MNRAS, 333, 481 

Beswick R., Pedlar A., Mundell C, Gallimore J., 2001, MNRAS, 325, 151 
Binney J., Tremaine S., 1987, Galactic Dynamics, Princeton University Press, 
Princeton 

Binney J., Tremaine S., 2008, Galactic Dynamics: Second Edition, Princeton 
University Press, Princeton 



16 



H. Engel et al.: NGC 6240: merger-induced star formation & gas dynamics 



Bland-Hawthorn J., Wilson A. S., Tully R. B., 1991, ApJ, 371, L19 
Boger G., Sternberg A., 2005, ApJ, 632, 302 

Bonacini Calia et al., 2006, in Advances in Adaptive Optics II, eds Ellerbroek B., 

Bonaccini Calia D., Proc. SPIE, 6272, 627207 
Bondi H., Hoyle E, 1944, MNRAS, 104, 273 
Bonnet H., et al., 2004, The ESO Messenger, 117, 17 
Bournaud E, Due P. -A., Emsellem E., 2008, MNRAS, 389, 8 
Braine J., Davoust E., Zhu M., Lisenfeld U., Motch C, Seaquist E., 2003, A&A, 

408, L13 

Brandl B. R., et al., 2009, ApJ, 699, 1982 
Calzetti D., et al., 2000, ApJ, 533, 682 
Capellari M., Copin Y., 2003, MNRAS, 342, 345 

Cappellari M., Neumayer N., Reunanen J., van der Werf P., de Zeeuw P.T., Rix 

H.-W., 2009, MNRAS, 394, 660 
Casoli E, Dupraz C, Combes E, 1992, A&A, 264, 55 
Chapman S. C, et al., 2003, Nature, 422, 695 
Chapman S. C, et al., 2005, ApJ, 622, 772 
Colbert D. T., et al., 1994, ApJ, 436, 89 
Condon J. J. & Yin Q. E, 1990, ApJ, 357, 97 
Condon J. J., 1992, ARA&A, 30, 575 
Conselice C. J., 2003, ApJS, 147, 1 

Cox T., Jonsson P., Primack J., Somerville R., 2006, MNRAS, 373, 1013 

Cox T., et al., 2008, MNRAS, 384, 386 

Cram L. E., 1998, ApJ, 506, 85 

Cresci G., et al., 2009, ApJ, 697, 115 

Dasyra K. M., et al, 2006, ApJ, 638, 745 

Dasyra K. M., et al., 2006, ApJ, 651, 835 

Davies R., Sternberg A., Lehnert M., Tacconi-Garman L., 2003, ApJ, 597, 907 
Davies R., Sternberg A., Lehnert M., Tacconi-Garman L., 2005, ApJ, 633, 105 
Davies R., et al., 2006, ApJ, 646, 754 
Davies R., 2007, MNRAS, 375, 1099 

Davies R., Mueller Sanchez E, Genzel R., Tacconi L., Hicks E., Friedrich S., 

Sternberg A., 2007, ApJ, 671, 1388 
Davies R., 2008, in The 2007 ESO Instrument Calibration Workshop, eds Kaufer 

A., Kerber E, (Springer- Verlag Berlin Heidelberg), p. 249 
Dejonghe H., 1986, PhR, 133, 217 
Di Matteo P., et al., 2007, A&A, 468, 61 
Di Matteo P., et al., 2008, A&A, 492, 31 
Downes D., Solomon P. M., 1998, ApJ, 507, 615 
Doyon R., et al., 1994, ApJ, 437, 23 

Drilling J., Landolt A., 2000, in Allen's Astrophysical Quantities, ed Cox A., (4th 

ed.; Berlin: Springer) 
Einstein A., 1915, Sitzungsberichte der Preussischen Akademie der 

Wissenschaften zu Berlin, 844 
Eisenhauer E, et al., 2003, in Instrument Design and Performance for 

Optical/Infrared Ground-based Telescopes, eds Masanori I., Moorwood A., 

Proc. SPIE, 4841, 1548 
Engel H., et al., 2010, ApJ, submitted 
Fakhouri O., Ma C.-R, 2008, MNRAS, 386, 577 
Falcon-Barroso J., et al., 2006, MNRAS, 369, 529 
Forster Schreiber N., 2000, AJ, 120, 2089 

Forster Schreiber N., Genzel R., Lutz D., Sternberg A., 2003, ApJ, 599, 193 

Forster Schreiber N., et al., 2009, ApJ, 706, 1346 

Gallimore J. E, Beswick R., 2004, AJ, 127, 239 

Genzel R., Lutz D., Tacconi L. J., 1998a, Nature, 395, 859 

Genzel R., et al., 1998b, ApJ, 498, 579 

Genzel R., Tacconi L. J., Rigopoulou D., Lutz D., Tecza M., 2001, ApJ, 563, 527 
Genzel R., et al., 2008, ApJ, 687, 59 

Gerssen J., van der Marel R. P., Axon D., Mihos J. C, Hernquist L., Barnes J. 

E., 2004, AJ, 127, 75 
Graham A. W., Worley, C. C, 2008, MNRAS, 388, 1708 
Greve T., Papadopoulos, Gao Y, Radford S., 2009, ApJ, 692, 1432 
Guilloteau S., et al., 1992, A&A, 262, 624 
Haring N. & Rix H.-W., 2004, ApJ, 604, 89 

Haring-Neumayer N., Cappellari M., Rix H.-W., Hartung M., Prieto A., 

Meisenheimer K., Lenzen R., 2006, ApJ, 643, 226 
Heckman T. M., Armus L., Miley G. K., 1987, AJ, 93, 276 
Heckman T. M., Armus L., Miley G. K., 1990, ApJS, 74, 833 
Hernquist L., 1990, ApJ, 356, 359 
Hernquist, L., 1993, ApJ, 409, 548 
Hibbard J. E., Mihos J. C, 1995, AJ, 110, 140 
Hopkins P. E, et al., 2006, ApJS, 163, 1 
Hopkins P. E, et al., 2008, ApJS, 175, 356 

Ikebe Y, Leighly K., Tanaka Y, Nakagawa T., Terashima Y, Komossa S., 2000, 
MNRAS, 316, 433 

Jogee S., 2004, in Alloin D., Johnson R., Lira P., eds, AGN physics on All Scales, 

Springer, Berlin, ch.6 
Johansson P. H., Naab T., Burkert A., 2009, ApJ, 690, 802 



Joseph R. D., Wright G. S., Wade R., 1984, Nature, 311, 132 
Karl S. J., et al., 2008, AN, 329, 1042K 
Karl S., et al., 2010, ApJ, 715, 88 
Kauffmann G., White S. D. M., 1993, MNRAS 261, 921 
Khochfar S., Burkert A., 2006, A&A, 445, 403 
Kleinmann S., Hall D., 1986, ApJS, 62, 501 
Knierman, K. A., et al, 2003, AJ, 126, 1227 
Komossa S., Schulz H., Greiner J., 1998, A&A, 334, 110 
Komossa S., Burwitz V., Hasinger G., Predehl P., Kaastra J., Ikebe Y, 2003, ApJ, 
582, L15 

Kormendy J., 1987, in Nearly Normal Galaxies, ed Faber S. M., (Springer- 

Verlag, New York), p. 163 
Kormendy J., Sanders D. B., 1992, ApJ, 390, 53 
Kormendy J., Illingworth G., 1982, ApJ, 256, 460 
Kormendy J. & Kennicutt, R. C, 2004, ARA&A, 42, 603 
Krajnovic D., et al., 2006, MNRAS, 366, 787 
Landini M., et al, 1984, A&A, 134, 284 
Lester D., Harvey P., Carr J., 1988, ApJ, 329, 641 
Lester D. F. & Gaffney N. I., 1994, ApJ, 431, 13 

Lutz D., Sturm E., Genzel R., Spoon H., Moorwood A., Netzer H., Sternberg A., 

2003, A&A, 409, 867 
Maiolino R., et al., 2007, A&A, 468, 979 
Mamon G. A., Boue G., 2010, MNRAS, 401, 2433 
Marconi A., Hunt L., 2003, ApJ, 589L, 21 
Max C. E., et al., 2005, ApJ, 621, 738 

Max C. E., Canalizo, G., de Vries W. H., 2007, Science, 316, 1877 
Mazure A., Capelato H. V., 2002, A&A, 383, 384 

Medling A., Ammons M., Max C, Davies R., Engel H., Canalizo G., 2010, ApJ, 
submitted 

Menendez-Delmestre K., et al., 2009, ApJ, 699, 667 

Mihos C, Hernquist L., 1996, ApJ, 464, 641 

Moriondo G., Giovanelli R., Haynes M., 1999, A&A, 346, 415 

Mueller Sanchez E, Davies R., Eisenhauer E, Tacconi L., Genzel R., Sternberg 

A., 2006, A&A, 454, 481 
Murray N., Quataert E., Thompson T. A., 2005, ApJ, 618, 569 
Naab T., Burkert A., 2003, ApJ, 597, 893 
Naab T., Jesseit R., Burkert A., 2006, MNRAS, 372, 839 
Netzer H., Lemze D., Kaspi S., George I., Turner T., Lutz D., Boiler T., 

Chelouche D., 2005, ApJ, 629, 739 
Netzer H., Trakhtenbrot B., 2007, ApJ, 564, 754 

Neumayer N., Cappellari M., Reunanen J., Rix H.-W., van der Werf P. P., de 

Zeeuw T., Davies R., 2007, ApJ, 671, 1329 
Ohyama Y, Yoshida M., Takata T., 2003, AJ, 126, 2291 
Origlia L., Moorwood A., Oliva E., 1993, A&A, 280, 536 
Pasquali A., de Grijs R., Gallagher J., 2003, MNRAS, 345, 161 
Peletier, R. E, et al., 2007, MNRAS, 379, 445 
Pollack L., Max C, Schneider G., 2007, ApJ, 660, 288 
Pope A., et al., 2005, MNRAS, 358, 149 

Rabien S., Davies R., Ott T., Li J., Abuter R., Kellner S., Neumann U, 2004, 
in Advancements in Adaptive Optics, eds Bonaccini Calia D., Ellerborek B., 
Ragazzoni R., Proc. SPIE, 5490, 981 

Robertson B., Bullock J., Cox T., Di Matteo T., Hernquist L., Springel V., 
Yoshida N., 2006, ApJ, 645, 986 

Rovilos E., et al., 2005, MNRAS, 359, 827 

Sanders D. B., Mirabel I. E, 1996, ARA&A, 34, 749 

Sanders D. B., et al., 1988a, ApJ, 325, 74 

Sanders D. B., et al., 1988b, ApJ, 328, L35 

Sanders D. B., 1992, in Relationships between active galactic nuclei and starburst 
galaxies, ed. A. V. Filippenko, ASP Conference Series (ASP: San Francisco), 
vol. 31, p. 303. 

Schmitt H. R., Bica E., Pastoriza M. G., 1996, MNRAS, 278, 965 
Shier, L. M., Fischer, J., 1998, ApJ, 497, 163 
Silge J., Gebhardt K., 2003, AJ, 125, 2809 
Silge J. D., et al., 2005, AJ, 130, 406 
Smail I., et al., 2002, MNRAS, 331, 495 

Soifer B. T., Neugebauer G., Houck J. R., 1987, ARA&A, 25, 187 
Solomon P. M., Downes D., Radford S. J. E., Barrett J. W., 1997, ApJ, 478, 144 
Springel V., Hernquist L., 2003, MNRAS, 339, 289 
Springel V., 2005, MNRAS, 364, 1105 

Springel V., Di Matteo T., Hernquist L., 2005, MNRAS, 361, 776 
Steffen A., et al. 2006, AJ, 131, 2826 
Sternberg A., 1998, ApJ, 506, 721 

Sternberg A., Hoffman T., Pauldrach A., 2003, ApJ, 599, 1333 

Sugai H., Malkan M., Ward M., Davies R., McLean I., 1997, ApJ, 481, 186 

Tacconi L. J., Genzel R., Tecza M., Gallimore J., Downes D., Scoville N., 1999, 

ApJ, 524, 732 
Tacconi L. J., et al., 2008, ApJ, 680, 246 
Tecza M., et al, 2000, ApJ, 537, 178 



17 



H. Engel et al.: NGC 6240: merger-induced star formation & gas dynamics 

Teyssier R., Chapon D., Bournaud, R, 2010, ApJ, 720, 149 

Thomas J., Saglia R., Bender R., Thomas D., Gebhardt K., Magorrian J., 

Richstone D., 2004, MNRAS, 353, 391 
Thornton K., et al., 1998, ApJ, 500, 95 

Tokunaga, A., 2000, in Allen's Astrophysical Quantities, ed Cox A., (4th ed.; 

Berlin: Springer) 
Toomre A., Toomre J., 1972, ApJ, 178, 623 

Toomre A., 1977, in Tinsley B. M., Larson R. B., eds, Evolution of Galaxies and 

Stellar Populations, Yale Univ. Obs., New Haven, p.401 
Tremaine S., et al., 2002, ApJ, 574, 740 
van der Marel R. P., van Dokkum P. G, 2007, ApJ, 668, 738 
van der Werf P. P., et al., 1993, ApJ, 405, 522 
Veilleux S., Kim D.-C, Sanders D. B., 2002, ApJS, 143, 315 
Vignati P., et al., 1999, A&A, 349, L57 
Wetzstein M., Naab T., Burkert A., 2007, MNRAS, 375, 805 
Whitmore B., McElroy D., Tonry J., 1985, ApJS, 59, 1 
Winge C, Riffel R. A, Storchi-Bergmann T., 2009, ApJS, 185, 186 
Wolf J., et al., 2010, MNRAS, 406, 1220 
Woo J.-H., Urry C. M., 2002, ApJ, 579, 530 

Yun M. S., Hibbard J. E., 2001, in Hibbard J. E., Rupen M., van Gorkom J. 
H., eds, Gas and Galaxy Evolution, ASP Conference Proceedings, Vol. 240, 
p.866 



18 



